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ABRSTRACT 
AM analysis was done of central 4rctic Ocean acoustic dates 
to determine the temporai and soatial characteristics of 


transient noise events. Digital ambient noise data from 
the FRAM IV experiment of April 1982 were searched for 
ambient noise transients using a detection proaram. The 
time series of the resulting detections were examined 
visually to categorize each detection as a transient. 
artifact or false alarm. The transient events were lacatec 
in Space using time delays between signal arrival at 
different hydroshones. The cross shape of the FRAM IV 
horizontal array permitted Iocation in both bearing anc 
range. The source strength of each event was calculated 
using a Simole dipole source model. Refraction and 
scatterring of the acoustic path in the Arctic Ocean was 
taken into account. 


ill 


TRe overail number af events detected, ana hence treir 
amnterarrivai times and spatial censity, were ail aftected 
Gy the backaround ambient noise level. The detection 
croaram used the same thresencid sicnali-tc-noises levei ts 
alli data tapes, sa when ambient noise levels were low more 
detections occurred. The mean interarrival time between 
events was 100 seconds. The interarrival time fit a J 
shaped gamma probability distribution. The number of 
events detected per area decreased with range from the 
array center. Half of the events occurred within 2006 
meters of the array. In this erea there were 9.2 events 
per square kilometer per hour. The event population showes 
mo predcminant anguiaer dependence. The strengths 
caiculated using the simple dipole model had a mean oF 439 
eM overall and 260 EN during Guiet times. Stronger events 
occurred during times with high ambient ncise levels. 
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CHAPTER i 
INTRODUCTION 
Thesis Motivation 

This thesis investigates the spatial distribution, 
strength and rate of occurrence of low frequency central 
Arctic noise events. During the last 30 years there has 
been a growing commercial, military and academic interest 
[eeeee eeectic region. 

In this relatively unexplored area there is increasing 
evidence of rich mineral and petroleum resources. Research 
into methods of locating these assets and constructing 
facilities to exploit them have received much attention. 
These racilities must be abie to withstand the harsh Arctic 
surroundings. The study af Arctic acoustics helps in. 
understanding the Arctic environment and climate. It has 
been shown that there 18s a direct correlation between i to 
2O Rez ambient sound pressure and environmental stresses and 
momentsliod. THe ability to use acoustic noise leveis as 
an environmental predictor would be a useful tool in the 
protection of commercial Arctic facilities. 

The Arctic ocean serves as a military arena for 
several submarine fleets. The underice environment makes 
detection difficult, increasing the stategic role of those 
fisets. Because af the sound velocity profile of the 
central Arctic there is a surface duct which channels sound 


for long distances. But, the underice profile scatters 





sound energy, effectively filtering out high 
freaquenciesCfiid. The result is that only low frequency 
Signals travel far in the Arctic, and therefore the low 
frequency range is the best for detecting adversary 
submarines. The eee tate of understanding the low 
requency ambient noise field becomes apparent. The actual 
central Arctic ambient noise level is at times much quieter 
chan the apenm ccean, But it contains unpredictable 
transient noise events which interfere with conventionai 
detection schemes. It has been hypothesized that the 
background ambient noise is the summation of these 
transients from throughout the Arctic basin£7]. Analysis 
of the spatial and temporal distribution of these 
transients is a logical next step in understanding low 
frequency noise, and improving our submarine detection 
capabiiity. 

TRe academic challenge of the Arctic lies in the 
sparseness of field data. The Arctic cannot be casuaiiy 
sampled. Even simple experiments require expensive 
expeditions. The harsh environment takes its toll on 
researchers and equipment, and reduces the amount of usable 
data. Hence, the study of the Arctic 185 like a jigsaw 


puzzle with few pieces oresent. The total picture remains 


t 


2 Stimulating mystery. 


This study analyzes data collected during the FRAM [IV 


ID 


“xperiment by Massachusetts Institute of Technology and 





Woods Hole OUceanocaraphic Institute personnel. The FRAM IV 
2ce camp was located in the Harents (Narmsen? Abyssal Flain 
at approximately 94° N by 15° E, as shown in Figure 1-1. 
The ice was S$ meter thick multi-year pack ice. The ice 
activity was lows there was no ice ridging or lead 
formation around the camp during the experiment. 

The FRAM IV ice camp was set up from 25 March toa Iii 
May if82. i Tnis study analyzes data taken between March 
27th and April 22nd. The weather was mild, with 
temperature ranged from —-25° to -4° C, and wind speed from 
(oven Ket Ss. 

The ambient noise was sensed with a large horizontal 
hydrophone array which consisted of two non-uniformly 
spaced line arrays, crossing at rignt angles. The data 


were digitally recorded on a multichannel system. 
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Figure i-it Location of the FRAM IV Arctic experiment 
Sememecee in the spring of 1982. Ciidl 





Thesis Contents 

My work in this area began long after the FRAM IV 
ambient noise data were collected. The first step was 
finding the events in the raw data. The data consisted of 
magnetic tapes each containing 20 minutes of digitized 
noise levels. A program was written which searched the 
ambient noise tapes for possible events. Chapter 2 
discrites the automated and manual techniques used ta 
accomplish this detection. 

These events were then located in space usina the 
difference in arrival time between hydrophones. This was 
also done with a computer program. The program plotted the 
arrival time delays against range to a trial location, did 
a least squares fit, and chose the location with the best 
Teens 15 covered im Chapter 3. 

The peak voltaces for each event were used with 
the dicole source modei to predict peak source strenagtn. 
Re background ambient noise strenath was also determined. 
These strength calculations are found in Chapter 4. 

In Chapter 3 the distribution of event interarrival 
time was determined. The event locations and strengths 


were analyzed, and a spatial density found. 


Chapter 6 summarizes the key results of this study. 





CHAPTER 2 
DETECTION OF NOISE EVENTS 

Data Coilection 

Twenty nine FRAM IV ambient noise tapes were searched 
in order to find a population of noise events for this 
study. The specific tapes were chosen from the possible 67 
in order to cover the entire range of days of the FRAM IV 
experiment. However, there were several days when no 
ambient noise tapes were recorded. To help fill these gaps 
five reverberation tapes were also searched. These tapes 
were recorded prior to the reverberation shot being fired, 
or they were recorded so late in the experiment (80 minutes 
after the shot) that reverberations were no loneer scresent. 

The FRAM IV experiemnt used the horizontai array oar 
omnidirecticnal hydrophones pictured in Figure 2-1. The 
Rydrophores were suspended from the ice into the water coo a 
Gepth of 9753 meters pelow the air/ice interface. ne two 
crossed iines of the array allowed the possibility cr 
localizing events in Space. Although 25 hydrophones are 
shown in Figure 2-1, only 24 at a time could be used toa 
record data. In most cases a few of the recording channels 
were used for other sensors (geophones or Rydrophones used 
ieee ver tical attrav). Most of the time 19 to 241 horizontal 
array Hhvdrcokone cata were recorded. 

The FRAM IV ambient noise tapes were rescorcead 


jest 


i 


Peaily. Figure 2-2 shows a scnematic of the system used 
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Figure 2-1 FRAM IV horizontal hydrophone array. [1] 





FRAM IV DATA ACQUISITION 
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Figure 2-2 Schematic of the recording system used for PRAM 
Po deea collection. (41 





tor coilecting the data. The input from each hydrophone 
was put through a gain ranging amplifier and a lcw pass &6 
Hz filter. Tt was sampied at 250 Hz and recorded on 26 
minute magnetic tapel4]. The 24 channel recorder had a 120 


qB dynamic ranqgeLliZ]. 


Svemeeetector Fregram 

TRe Svent Getector program was written to uake digitel 
data from a FRAM IV tape and determine where in that tape 
noise events occurred. The program was originally written 
to take data cirectly from a tape drive, But subseauentiy 
modified to take the data from a file. The framread 
Program, with a -head switch is used to read the tape into 
tne file. This wiil eliminate any headers and then reac 
the digital data straight into a file. A FRAM II tape mav 
b= read into 3 similar file using framread and the sweitcnes 
-head and -fram2. The framread procram was written by 
S. Duckworth, and is available to the Arctic Acoustics 
Frogram at MIT. 

The event detection programs source codes, ficw chart, 
and a short users guide are found in Appendix A. The event 


detection program which reads from a file is called 


hdetect. The detection proaqram reads from a tape crive. 


We 


Both procrams are written in the c grogramming language far 


a UNIA operating system. 


De 


Tre esvemt detection program follows the Bbiock diagram 





— 


Soe iaqure 2-3. The program initialization portion defines 
variables and constants, zeros flaqs, and requests user 
inAput such as tape number, date, time and 

channels (hydrophones) to be used, as well as, the name of 
input and output files. Arter this information is 
requested from the user, the program no longer requires 


Se eent.s OM. 


Re event detection ocrogram then reads in a file af 


data, filters the data, squares each data point, and takes 
the square root. The filter was a Parks-McClellan digital 
“=O to 80 Hz bandpass filter. Its frequency response 15 
Shown in Figure 2-4. The range of this filter was chosen 
to avoid the Nyquist frequency (125 Hz) and hydraphone 
cabie strum (1-20 Hz), and to be compatible with the’ anaicag 
S30 Hz low pass filter the data went through before being 
recorded. The data were squared and then rooted to ensure 
DOSitive oeak vaiues for ali Gata points. 

tae 2ext portion of the program used 4 thresnoisa 
Getection scheme to check each channel for possible noise 
events. For a particular channel a short average of the 
four most recent data points was compared to a lonq average 
Of 64 recent data points. If the ratio of short averace toa 
long average was over a certain value that channel would be 
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Lagged for 4a possible event. The time of the flaaq and the 


ul 


value of the shert average were also recorded. ALL other 


channels averages were taken similarly. If at least SOX of 
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Figure 2-2 Block diagram of major modulues of the 
detection proaram. 
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Figure 2-4 Frequency response of the Farks-McClellan 
digital bandpacs filter used in detection. 
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the channels were flaaqged, an event would be declared. 
Then the program would shift to the next time increment of 
data and the process would be repeated. 

A more detailed diagram of the event detection moduie 
is seen 1n Figure 2-S. There are four Bee reset 
flag, set flag, new event, and deactivate old event 
modules. The reset and set flag modules deal with the 


Channel flags which trip when a particular channel 


experiences a larage signal-to-noise ratio (i.e. the RATIO 
of short average to long average exceeds a certain level). 
The new event and deactivate old event modules deal with an 
active event matrix which identifies active events, and 
stores channel flag time and amplitude for each cgeclared 
event. 

The reset flag module resets the channel flaa if it 
has been more than ©.3 seconds since the channel tripped. 
Spurious peaks on a channel might flag a channel 
Orematurely. This reset module prevents a number of 
channels with spurious peaks over a long time period from 
being falsely deciared an event. The value of this RESET 
DELAY was determined by examining known events and noting 
that about half the channels tripped within a 0.32 second 
period. 

The set fiag module determines 17 the short averace ta 
long average RATIO has been exceeded and, if it nas, the 


module 1) checks to see if the channel flaad is already 
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Figure 2-S Diagram of event detection module flow and 
decision making. 
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BereGed, =<) checks ta see if this detection is part of a 
recentiy declarec event, and 3) adds its detection 
information to the active event or the channel flag, or 
discards tne information depending on the circumstances. 
The short average length ct 4 (= 0.016 secd was cnosen so 


that 1t matched the length of the signal (9.02 sec). This 


provides the maximum signal level since this is long Ssnough 


ccs Gen ail of the signal ane short enough mot =O average it 
with the tower surrounding backqround noise level. the 


tong averace length af 64 was chosen because a Llenath ratia 
Of 15 to 1 Rad been suggested by KFellyC€91] for the Large 
Aperture Seizmic Array (LASA). The detection and 
localization schemes used bv this Llarae horizontai array 


if 
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Were wWirectiy applicable to tne FRAM {Vv Bydrophone array 


Gata The choice of detection RATIO was done through a 
Series of tests. anc the selection wes made By Halencine 
Getection rate anc false aiarm rate. 


Wren an event 1s ceclareaG the intormation in tne 
Channel flags is transferred to the active event matrix, 
and tne channel flags are cleared. So if a channel has a 
detection and its channel flag is not already tripped. the 
set flag module must first see if the detection beiongs to 
a recently deciared event. If the active event already nas 
shat channel flaccgeac. the information is replaced ov the 
mew detection aniy if their times differ by iess than O.02 


seconds and the new detection amolitude is greater. This 





means that the existing data can oniy be replaced by a 
aetecti1o0n af the same signal Raving a higher peak value. 

If the active event does not have that channel flagged, the 
Getection intormation is entered in the active event 
matrix, and the channel flag is not tripped. The oldest 


active events are checked first, and the detecticn 


information entered in only one event. If all of the 
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as been more than 0.02 sec since the most recent inout, 
the detection is considered potentially part af an 
unceclared event, and its channel flag is tripped. 

If a detection is made on a channel and its channei 
flag is already trioped, the information in the channel 
flag will be replaced with the new cetectian oniy if the 
new detection amolitude is greater. 

TRe new event mocule checks to see if at ieast SO" at 
the channel flaqse have besn tripped, and if so. deciares a 
new event. The SQ%Z mutual occurrence criterion was used if 
the LASA program with good results€S)]. All the channel 
flag information is transferred to the active event matrix, 
and the channel flags are reset. 


The deactivate old event module was used to remove 


events which were pest. A set EVENT DELAY time atter an 
e.ent is deciared. 14 165 written to tne outout Ffiie ana 


0 


rased trom the active event matrix. This prevents 


wious peaks from being added to a event ionag past. TRe 


Wn 
OQ 
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EVENT DELAY time of 9.5 second was chosen because 
inspection of the known events revealed that most channeis 
tripped within a 0.8 second period. The RESET DELAY (0.03 
sec) plus the EVENT DELAY result in this 9.8 second look 
for each event. 

The output of the event detection program is an output 
file which contains the time each event was declared, which 
channels were +lagged and the time delay and peak amplitude 
for each channel. The time delays are relative to the 
earliest channel flag time, so one of the channels always 
RhaS a zero time delay. 

The first version of the detection program was written 
tG take a new short average and long average at every data 
peint ‘every 0.004 seconds). TRis program took 4 toa 8 
hours to search a 29 minute data tape. A concession toa 
speed was made and the program changed te computs averages 
at every fourth sample point (every 0.015 seconds). This 


reduced th acy of the time delays anc the abliity a7 
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fy 
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the program to pick up events. The RATIO had to be lowered 
im order to get the same detections which were obtained 
previously. 

Studies to find the best signal-to-noise RATIO were 
conducted several times. Development of the LASA detection 
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system nad reveasied that a 7 AB signal-to-noise ratia was 
oe 


needed for 75% detection(al. This equates ta a 3 ta l 


ratio of siqnal power to noise power. Since [ was working 





with pressure vice power I used 2.2 as my starting RATIO. 
This RATIO engulfed me in false alarms. A quick study was 
done around the 2.4 level. The first 10 minutes of tape 
4013 were run at RATIOS of 2.4, 2.45 and 2.48. This tape 
had been visually examined in detail previously, so the 
events were known. The results are shown in Table 2-1. 


Slso shown in Table 2-1 are the results of a second study, 


eone after the program had been changed to average ies 
often. 
Table 2-1 Determination of the Best RATIG 

Average taken at every data point 

RATIO Fee Stan wre te false alarm rate 

2.4 P47 fe ctite 

Bato ats ee 

2.48 TE ee 
Average taken at every fourth data point 

RATIO detection rate false alarm rate 

aS Tor 457 

a oo rales 7 

a 59% Boi 

fa a7 aoe 
Getection rate = # event detections / # of known events 


FA rate = # non-event detections / total detections 





Notice how the detection rate has decreased and the 


false alarm rats increased as a result of only averaging at 


ae 
ye 


every fourth data point. Averaging less frequently means 
there iS a smaller probability that the deta points to Fe 


averaged will ail lie near the peak amplitude of the 


si 


if} 


Mal. The signal level is generally lower than that 


detected when averaging every data point. and a iower 


Signai-to-neise RATIO must be usec to detect thse same 


pols 
rly 


events. But when the signal-to-noise RATIO is lowered the 
false alarm rate increases. 

The RATIO of 2.28 was settled on. This is a 
compromise which gives a detection rate which finds most 
high and medium strenath events. and which has a tolieratis 
faise alarm rate. Because the detection rate is less than 
100% (714) there were events present which coulc be seen 


Yisusliy, But were not sicked uc by the detection orscaranm. 


The RATIG could have been adjusted to detect ali events 
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seen visually, but at the cost of a multitude cf Fail 
alarms. The RATIO was kept at 2.28 and used for the 
Getection of all data tapes. 

The final version of the Adetect program read digital 
data from a framread file, detected possible events usina 


the less fr 
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quent averacing scheme, and supplied the event 
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ead 


avs and amplitudes ta an output 
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Time anc chann 


rhe 


file, Gnce the RATIO had been satisfactorily set the 


program was used to search the FRAM IV ambient noiss tapes 





ror events, and no further program development was done. 
Weaknesses in the program were subseauently discovered, But 
have not been corrected. 

The biggest problem is the accuracy of the event time 
(the time when 30% of the Rennie have been flagged). An 
event time reported by the event detection program will not 
exactly match that found by plotting the time series. The 
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ths 


re usually within 3 seconds of each other. out have 
Been off By as much as 13 seconds in one caste. TRe time 
difference between the two methods is greater at the end af 
a tape, and is iikely toa occur after a particularly strang 
event has taken place (though there were times when time 


discrepancies develosed without strong events present. anc 


¢ 


aiso many strong events existed which dic not induce 


discrepancies). Typically, there might be ro time 


ak + 


iti 


rence at the first part of the tape, then after a 
strong event a three second discrepancy would se seen ane 
182 S would Ge consistent until the end of the tase. 
secause the errors did not appear randomly throughout the 
tape, and because they developed impulsively. I believe 
that the problem lies in the time counter of the event 


detection program becoming offset from the time of the raw 


data, perhaps because of short records in the raw data. 
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v tection goroqram was mot written £0 


time keener in the event of a short record. Correcting 


this may eliminate the time discrepancy oroblem. 





It has already been mentioned that the accuracy of the 
time delays deteriorated when the program was changed toa 
run more quickly. This also made the detection of the 
event peak amplitude less likely. As a result, in order to 
get accurate locations and source strengths, both time 
delays and peak voltages had to be taken manually from time 
series plots of each event. 

The cther major problem of the event detector program 
1s that it does not discriminate between an Arctic noise 
transient and an artifact such as an air gun blast ora 
reverberation shot. Short, strong signals are reported as 
possible events. Adding this discrimination to the program 


is the next step in improving its usefulness. 


Visual Confirmation 

Visual confirmation was required for all possible 
noise events in order to eliminate artifacts, false alarms 
and multiple detections of the same event. In addition, in 
a few cases visual confirmation revealed two events where 
there had only been one detection. 

The event detection program was designed to preclude 
the need for plotting a time series of each event. fhe 
output of the program contains time delay and voltage 
amplitude information which can be used directly in the 
location program. However, because of the decreased 


accuracy of the time delay and amplitude information, and 





because of the event time discrepancy mentioned previously, 
it was necessary to plot the time series of-.each event. 

The first step of the visual confirmation is to review 
the tape log for any artifacts that may have occurred 
during the recording. The times are noted, and these are 
compared to the event times given by the detection program. 
Then a time series of the artifact was plotted to determine 
which detections were associated with it. In general, an 
artifact such as an air gun blast did not affect detections 
for over 20 seconds. 

The visual confirmation portion of the procedure 
evolved from a very limited look only at events which could 
not be located with the detection program generated time 
delays, to a three step plotting procedure for each event. 
During the early period of this work the Adetect program 
Output was used directly as the input to the location 
program. The location program used the time delays to 
determine the event’s location in space. Those events 
which could not be located needed a closer look, and so 
their time series were plotted. The plots were made of a 2 
second period including the event time given by the hdetect 
program. Often there was no apparent event in this time 
series plot, and the detection was declared a false alarm. 
When an event did plot, manual time delays were taken and 
used to locate the event. These manual delays located 


these events with better accuracy then the hdetect 





Generated time delays. Tt soon Gecame apparent that the 


cest answers would be obtained by taking manual time cGelays 


Se 


of all events. Trying to localize the events with the 


program generated time delays was drapped from the 
procedure, and the first step after getting the Adetect 
program output Ssecame doing a 2 second time series oloat of 
each event. 
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the discovery of the event time discrepancy was made. 

After plotting a dozen events right at the time shown by 
the program, the final two dozen event of tape 4009 all 
appeared to be false alarms. The quality of the tape was 
good (‘low background noise), so this seemed hichly 
Suspicious. A Droader search of the time around each event 


showed that the final two dozen events were not false 


alarms, but were events with times 3 seconds different than 
those indicated by the program, so that none of those 
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events had shown up on the 2 second time series pl: 


Ul 


method of visually confirming events was cnanged so that a 
waterfall plot was made around the time of each event. 

This showed the exact time of the event, and helped discern 
the pattern of time discrepancy between time series and the 


Adetect oeroegram. Once the pattern was found events were 
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confirmed events were then plotted with 2 seconde time 
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jTRe finai stes of visual confirmation was simply 
plotting and replotting the events to the proper aain sc 
that the higher voltage amplitudes would not be cut off. 
TRese peak voltages were manually taken from the time 
series plots and used to determine source strength. 
The final method used for visual confirmation was: 

1) Check tape log for artifacts, and eliminate 

those from further analysis. 


2) 


= Fiot a waterfall time series around each 
possible event, separate real events fram false alarms, and 
find true event times. 
ae? Flot = second time series of each real event. 
adjusting gain to keep from ciipsinag higher voitages. 
Using this technique certainly reduced tne false aiarm 


rate. A Breakdown of the detection statistics for tapes 


trat hed been examined BY Both methods is found in 


Table 2-2 Detection Statistics for Two Visual 
Confirmation Methods 


Artifacts Events False Alarms 
Orignai Method 16.94 Ba eh de 45.67% 
wif 2S sec plots 
Lliltimate Methoa L6n 2% 65.424 14.4% 
w/ Waterfall siot 
oe See Fi ok 


(Fercentagas of getections classed in each cateaory) 
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The human interpreter was a necessary tool in this 
scheme. There was not necessarily a one-to-one 
correspondence between events and detections. There were 
cases where a strong event would cause multiple detections, 
and cases where two events occurred at the same time and 
caused only a single detection. In some cases a series of 
detections seemed to be an event and an echo, or perhaps 
Straf. This would be counted as a single event. 

The method of determining whether a detection was a 
false alarm or a weak event was sometimes difficult. In 
general, if the detection program indicated a possible 
event, "something" could be seen on the waterfall plot. 
The detection was dismissed as a false alarm if no’ pattern 
for taking time delays could be seen. (Because of the 
Shape of the hydrophone array there were consistent 
patterns of time delays depending upon the direction to the 
event.) Presumably the false alarm rate depends upon the 
training and attention of the human interpreter. 

Manual time delays were taken from an arbitrary 
reference to the crossing of the largest peak to peak 
amplitudes, as shown in Figure 2-6. For most events this 
was Clear, but for weak or complex events some intuition 
was needed. 

Voltage amplitudes were taken as the maximum peak 
voltage in the event signal. All were taken as magnitudes 


regardless of sign. 








Figure 2-6 Sketch showing point of measurement for event 
time delays. The event is timed at its zero crossing 
between the largest pair of positive and negative peaks. 
(The measurement/analysis system has a polarity of negative 
voltage for positive pressure.) 





Ail types of noise event signatures previously 
observed by Dverl7] were seen in the ambient noise tapes I! 
evaluated. The majority of events were pops and extendec 
pops. There were aiso a few whines and strat events. 


While signature types were noted in general, the signature 


type of each individual event was not recorded. 





CHAPTER 3 
LOCATION OF NOISE EVENTS 
Event Location Frogram 
The program used for localization was based on the 
program FRUAK by Feter Stein Cit]. This program places the 
event at different trial locations and computes the slant 


range to each Aydrophone. Figure 2-1 shows the coordinate 


un 


ystem used far these Caiculations. TRese are opiatted 
against the time delays and a least square fit is done to 


determine slope as shown in Figure 3-2. 
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y intercept = 


The standard deviation of the time delays from tre sicre 


Line is figured. 
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The location having the lowest standard deviation is the 
location of the event. The inverse of the siope is the 


Group speed of the signal. The y intercept of the eplot is 


My 


dded to thre reference time of the manual time delays 


rr 


Rex 
get the time the event actually occurred (as opposed to 


when it reached the hydrophone array). 








Figure 3-1 Coordinate system used for calculation of event 
Ree ct2 OM . 
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Figure 3-2 Least squares fit of time delays and test 
location slant range. 





I have assumed that the strongest peak pressure sensed 
at the hydrophone is due to a waterborne acoustic 
propagation path from a source located in the ice sheet. 
The signal enters the water near the source and propagates 
directly toward the hydrophone. 

IT have assumed that the signal does not bounce off the 
Ocean bottom or the ice canopy before reaching the 
hydrophone. Faths bouncing off the bottom would produce 
signals with much lower energy than the direct path signal, 
and can be ignored. Signals bouncing off the ice canopy 
are too energetic to ignore but, as I show subsequently, 
they do not affect the time delay computations significantly. 

The location program is based on arrival times being 
related to slant range, R , and does not take the upward 
refraction of the acoustic path into account. The impact 
this has on the results is discussed in the next section. 

The location program takes as input a file of time 
delays and voltage amplitudes, and outputs a file 
containing the best event location, sound speed, and 
standard deviation. It also computes source strength based 
on the voltage amplitude inputs, the event locaticn anda 
spherical spreading loss. This feature was originally 
included so that the source strength could be computed 
directly from the event setection pregram outputs. Since 
the peak volv.ages recorded by the detection program are not 


as accurece as those done by hand, and since the 





transmission loss does not follow simple spherical 
spreading, these computed source strengths were not used 
for any part of this study. 

The FeUAK program set up a grid of points around a 
specified center position. The @rid consisted Gf a soint 
every 100 meters from -—-5000 to +5000 meters in both the x 
ana y Girections. This resulted in 1900 x 100 test 


bocations. When the best test location was found the 


Ul 


interval spacing was reduced to every 109 meters, and 
another 10,000 test locations were generated using the best 
location of the first round as the new center. The orocess 
was repeated with a i meter interval to get the final 
answer. The scheme evaluated a total of 30,000 test 
Locations, covered a range out to 30090 meters. and cack 
about 2O minutes to run. 


I moted that 4&4 sianificant number of events found with 


FaUAK were at the range limit of S900 meters. The pragram 
£0Ocation was written to searcr a larger area taster. Toes 
Fineness of the grid was decreased to 290 x 20 vice 100 x 


100. A 1000 meter interval was added to enable the orcaram 
to search out to 19,9000 meters. This reduced the total 
number of test locations to 1600 (20 x 20 ~4), and the time 
to one minute. lecation gave answers which were very 
censistent with Fa@UYAK, except in one particular situation. 
The wider grid size led to one problem. The lLocatian 


program sometimes found the lowest standard deviation for a 





point in the quadrant directly opposite the true location. 
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his was sugcgested by the sound speed being reported as 


approximately -1440 m/sec, as illustrated in Fiqure 3-3. 
This problem was solved by modifying location to make the 
program fTineloecate. This program used the grid size and 
spacing of FAYAK, and centered the search so that the user 
could designate which of the guadrants would be searched. 
A casuai look at the manual times delays of an eSvent easily 
reveals the appropriate quadrant. This program works well. 
but is as slow as the original F@tUAK. It was used rarely. 
As with the original FRUAK, I began to notice that 
some events were located at the range limit of the lIecation 
program. Mivosied te the modification of the location 
program to rorm the program fariocate. This proacram uses 
the flocatren grid size and fineness, But aliows the center 
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coint to Ge any of the far corners of the cricinal ilaoecat 


Junt 
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Sr 1d or at the limit range at each of the cardiral points. 
—_ 3 = 2 


Tris i 


uf 


ehown in Figure 2-4. This ailowed events to fe 
located out to 20,000 meters. 

The lIecation program source code and a brief user's 
manual are found in Appendix &H. This program was written 


in the c programming language for the UNIX operating 


system. This program was developed to the point of 


mess, and chen used to lincate events. Mo further 


Uy 


program deveicoment was done fexcept the very minor changes 


ts produce finelecate and farlocate), so there are surely 
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Figure 3-3 Least squares fit to a false location 189° 
away. 
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Figure 3-4 Areas covered by location. farlecate and 
finelocate. 








improvements to be made. 

The location program is quite interactive. One 
hydrophone with a bad time delay can change the slope and 
lecation a great deal. An event is located by eliminating 
bad time delays and checking the sound speed and standard 
deviation of the location. In some cases no hydrophones 
needed to be removed, but in most cases at least one 
hydrophone was removed before an event was considered 
located. The sound speed was the major indicator of 
whether an event had been located. If the sound speed was 
between 1390 and 1500 m/s the event was considered ee ceae 
Of course an attempt was made to get close to 1440 m/sec. 
This had to be balanced with reducing: the standard 
deviation. A standard see on below 0.01 seconds was 
considered good. 

A&A table summarizing all of the events and their 
location parameters 1s found in Appendix C. The standard 
deviations ‘sigma) are given 1n two sets of units. The 
first is the sigma calculated by the location program, and 
1t 1S in seconds. The second sigma is a translation oft 
that standard deviation to meters using the sound speed 
calculated for each particular event. The standard 
deviations ranged from 90.0010 to 0.0327 sec, with 90.0077 
sec being the average. The significance af this standard 
deviation will be discussed in the next section. 


In some cases just removing suspect time delays did 





not lead to a localization. A reexamination of the event 
time series was done to see if any of the manual time 
delays was incorrect. Often a reexamination of the time 
series produced a change of i to 4 of the time delays. 
These corrected values plus values from the other channels 
would then be used to locate the event. About 15% of the 
events required reexamination. Most of those were 
subsequentiy located. 

Despite the above efforts, there were a few events 
that could not be located within the 1380 to 1500 m/sec 
sound speed limits. These events may be from propagation 
paths other than the assumed direct acoustic path. Events 
arriving primarily through the ice longitudinal wave or 
the ice flexural wave would Rave phase speeds above and 
below my sound speed limits. These non-locatable evvents 
ar2 indicated in the event location summary of Appendix C, 
and they were not used for any analysis which required 
accurate Location. 

Figure 3-5 shows the position of the events located 
within a 2 km square centered on the array origin. Figure 


5-6 shows the position of all events located. 





2° 
co wo — 


S 999900 
—_ OO Ga > HM «Ii 


(meters) 
(Thousands) 
Este to. ei el 
ooo go 9 9 2 2S © 
© CO = oo echo 


{ 
el 





= -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 ] 
(Thousands) 
(meters) 


Figure 3-S Noise events located Miehim@ a © en SQuere 
surrounding the array origin. 
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etftects of Retraction on Location 

A sound speed profile was used to get the refractive 
paths for various ranges. This was simplified by the fact 
that all aee Ahydrophones were at a depth of 93 meters. 
Assuming that enly the "direct" path is involved means that 
each horizontal range has to have a unique launch angle in 
order to reach the Rydrephone at its specific depth. Rays 


were iaunch 


ful 


G into the layers of the sound velocity orofiie 
and the horizontal range to the Rydrophone was caicuiated. 
The time required to travel the refractive path can be 
calculated and compared to that of the slant range. This 
time error can then be related to the error of the location 
procram. 

Figure 3-7 shows the tinearized sound velocity profile 
that was used. It i185 based on the sound velocity orofiie 
reported for the eastern Arctic ocean By ChenfldI. Figure 
3-83 helps to illustrate the scheme used to calculate tre 


<= 


ray paths. Equations 3.4, 2.3 and 2.6 were used to 


calculate angles, ranges, depths, and propagation time. 
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Figure 3-7 Sound velocity profile used in predicting 
refractive paths. 
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Figure 3-8 Scheme used to compute refractive propagation 
paths. Source is located at z = Om and hydrophone at z = 


93m. Sound speed gradients chance at 80m, 254m and 362m. 





Lsing the known depths (2) and estimated speed 
gradients fa) of Figure 3-7, and choosing a particular 
launch angle (SQ) > all ar the subsequent angies of 
intersection of the layer interfaces (4,, 8=, Gq) can be 
found from Equation 2-4. The angle which intercepts the 
hydrophone at its depth (8,5) can also be found. With the 


angles known. the horizontal range and propagation time for 


pa 
rt 
LL 


each - can b&b 


{Lh 


ay determined with Equations 3-2 and 3-4. 
These are combined to get the total HRorizontai range and 
propagation time. Tt should be noted that there are two 
ranges and times for each launch angle. The first path is 
that which intercepts the hydrophone on the way down, while 
the other intercepts the hydrophone as it 1s refracted ‘back 
toward the surface. The maximum depth reached by the 
oragpagation path was also found, and those paths that went 
Sees meters. resulting im a range greater than 30600 
meters were not reported. A tabular summary was made ar 
launch angle, morizontai range from the hydrophone, maximum 
depth and propagation time, and this may be found in 
Appendix D, along with more detalled tables listing &,_,4; 
scat and ty _4- 

Rays connecting source and hydrophone with one cr more 
Bounces from the ice were not considered here. Tre eftfect 
of thos] rays will be taren into account in Chagoter 4. 

The refractive Gropagation time was iess than the 


Slant range propagation time because the refractec path 





travels through faster water. The slant range propagation 
time was calculated by dividing the slant range by 1428.48 
m/sec, the average sound speed between 0 and 93 meters 
depth. The time difference between the slant range path 
and the refracted path are shown in Figure 3-9 as a 
function of horizontal range. 


This time difference is greater than the average 
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o Toe ee eae ae nate deo ge ; 
G on of the iecation program aoniy after 


pede 
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Yianz 


CE 
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teangar 


ui 


1 


13,9000 meters, and the time G27 fer enee at 20,000 metercmis 
only about twice that average. The standard deviation does 
not reflect the time difference due to refraction because 
all of the hydrophone time. delays are adjusted in the same 


manner and direction. Figure 2-10 shows that siaqma doe 


li 


CL 


not grow with Horizontal range. Refraction effects Go not 
influence the standard deviation greatly. Closer than 


12,9090 meters the rance errcr caused by other factors mask 


WN 


any error fram ignoring retraction. 


etter point of focus for examining the 


tL 
Cr 


There is 
effect of refraction, and that is the change 1n time delay, 
not the change in the propagation time itself. A point was 
chosen at approximately 3S km from the origin of the 
hyaGrophone array, and another chosen at approximateiy 6 km. 


The slant range propagation times and refractive 


j 
nh 
labs 
(0 


ereoagation fines were calculated for each point. TRe tim 
Gelay between these two points was 9.6857 sec for the 


refractive path and G.6920 sec for the slant range path. 
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Figure 3-9 Time difference between the slant range 


acoustic path and the refractive path as a function of 
horizontal rance. 
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Figure 2-10 Scatter plot of standard deviation and range 


for all located noise events. 








The difference was 0.0075 (about 1%), or a deviation of 
QO.Q0Q057 sec for each hydrophone. Two points at 
approximately 18 and 20 km were also evaluated. The 
difference between their time delays was 0.01982 sec (about 
Q.14), or Q.0091 sec per hydrophone. These numbers are the 
same order as the total error of the location program. 

The location program may compensate for some of this 
error by raising the sound speed. If just the points above 
were used, the sound speed would go from 1438.5 to 1453.8 
at 3 km and to 1458.46 at 20 km. With 24 time delays being 
used in the location program the effect may not be as 
great. 

The main source of error in the location proqram is 
the quality of the manual time delays. When the signal-to- 
noise ratio was low, picking the correct peak was often 
difficult. The standard deviation wili reflect the 
judgement of the person picking off the time delays. The 
time delays were only measured to the closest 90.003 sec. 

Tt 15 interesting to note that 9.006 seconds equates to the 
width of a pencil tip on the time series plot scale. 

The final question to be answered is "How do the 
standard deviation and refraction errors equate to the 
range and bearing accuracy of the lecation program." Two 
nypothetical noise events were investigated, one at S000 
meters (2845, 4126) and the other at 20,000 meters 


Saseoso, 1S896) The time delays for slant range propagation 





and refractive propagation were calculated. The location 
program was run for each set of time delays, and for each 
set partially contaminated with 0.016 sec errors. (Zero, 
+0.916 and —0.016 were each added to one-third of the time 


delays.) The results are summarized in Table 3-1. 


Table 3-1 Results of location Program Accuracy Test 


R /\R O Ag Oo; c 


(m) Cm) (deg) (deq) (sec) (m/s) 

slant range 

=o km 4725 See 25.6 +O, 1 0.0003 1440 

20 km 18477 —-2143 DS iae/ eet Oar Q. 0002 1437 
refraction 

2 km 4724 ee Sele Oat 0.9003 14354 

20 km 18477, -2145 113.7 +0.1 0.0002 1458 
slant range w/ 

Sain 4191 =829 Stee et oS 0.0127 1453 

20 km 20526 —94 1 oS tO, S 0.9130 1443 
refraction w/ 

9 km 4209 et Soca} =tOns 0.90130 1446 

20 km 20526 —9G4 et eee Oo 0.90130 14623 


The refraction contaminated by errors case is closest 
to what was input into the Jocation program for the field 
events. This table gives an estimate of the accuracy ot 
the location program as 800 m at 3S km, and 2000 m at 20 km. 


The bearing accuracy 1s excellent. 





CHAPTER 4 
STRENGTH OF NOISE EVENTS 

AcCoUStTIC Saurce Model 

The dipole 1s considered a possible source model. 
Peak values for the source parameter of force, F , are 
used. 

The acoustic pressure due to a non-convecting compact 
Gipole source, in a nonrefracting infinite medium, isi: 
sin 8 fe Ghee 


ene) ; (4-1) 
477R fe Gt R | 


where RK = slant range. 


Figure 4-1 shows the orientation of the presumed disole. 
The angle aga is the launch angle from the horizontai osiane 
down into the water. 

Assuming that F may be exoressed as a harmonic. 
ld F/Ot| =(QJF = 27(fF . the pressure may then be expressed 


aS 


sin Sale 7(tr © F ; 
 ————— a | ee (4-2) 
47, R c R 
Solving for F gives 
ae 47( R“c 
P= a eo 
ne Smet ete | 


In the far field the 27(fFR term dominates the sum in 
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Fiqure 4-1 Assumed dipole source orientation and 


definition of launch angle, 9 . 





the denominator. For the lowest frequency considered in 


Ld 


2 — 


mewsetey (20 He), ¢ (1440 m/sec) is 10% of BUFR at 11 
meters and oniv 1.1% of 297 FR at 1000 meters. When the c 


in the cGenominator is nmeaqiected the force can be written 


aS. 


=p ok 


— (4-4) 
oO sin @ © 





where ) is the wavelength. 


The peak pressure, SR COUR em ead EO Ene (Peat ae otgae 


Pe - This definition of force was used as the parameter 


for dipole strength. Event signatures that were recorded 


trom a source within 200 m of a hydrophone were not usec t 


Calculate dipoie strength, F, : from peak pressure, p 


aoe 
For this medei the peak acoustic pressure must ba 
+ound. Tre Avydroohene sensitivity ot -is? dB re i voit pear 
Z [Les Was used to convert voltage to pressureli7s. 
oe 
Pvolte=- 67 N/m = 89 Fa . (a) 


The dicole strength formula requires wavelenath. IN 5 
Frequency was taken from the time series plots for each 


event vila axis crossing rate, and I was decermined by 


dividin c (1440 m/sec) By the frequency. 


3 —z = Fie cms. ae } a £ 8: “—< 
urch ancq.e needed tor the dipole model can 5S Founs 


tI 
i 


' 
* 
come *. 


as in Chapter 2 ty essuming a sound velocity proflis 
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computing the r 


Ta) 


Pave awit. There 15 a unique taunch 





angle for each horizontal range when the path is purely 
refractive, But a range of launch anaies when surface 
reflecticn paths are included. 

The final parameter in the dipole strength formula is 
Siant range. The answers obtained using slant range are 
the strenaths Based on spherical spreading in a non- 


absorptive medium, equation 4-4. Recause the spherical 


joo 
dU! 
lt 
Fi 
=r) 
y 
fl 


ing assumption is poor one, ‘refractive and sur 
reflective propagation paths are caused by the Arctic sound 
velocity profile), equation 4-4 must be modified. The 
effect of refraction on spreading loss will bs discussed in 
the next section. 

Volumetric absoretion was found By using the 
absorption formulas of Dyer(SJ. Assuming a GH of G.2, a 


Salimity of S35.5 “foo, a temperature of 0° C and a pressure 


of 40 atmospheres. [I calculated the total valumetric 


— 


ame Te 
we «9 


BH@seGeeciom to be 1.2 x 10 dqdBR/km tor an 80 Hz Siqnal. Far 
my maximum Horizontal range of 20 km, the absorption wouiad 
be 0.9026 dR. This 1S not significant, and I therefore did 


mot include a volumetric absorption correction in the 


strength calculations. 


Effects of Fefraction on Transmission Loss 


un 


cal spreading 1c6S in @& nOonretfractingag medium is 


da 


oher 


pels 
pea 
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rr 


ustrated in Figure 4-2C0141. The sound pressure 1s 


sressumed to sporead radially. The sound pressure squared 
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Figure 4-2 ‘Spherical spreading. (Sigal Umpiteibi4 |) 





is proportianai to intensity, and intensity is the power 


Cer unit area. Since the power from a source is constant 


The intensity at the reference range of 1 meter can be 


related to other intensities by 


Teo 47T I of 


vv 2 wr 


s+ (4-7) 
R > = " 
ayy R= Bee 





Since [ = p*/ Oc > this can be expressed in terms af 


transmission loss, H . 


tJ 


H = -10 log ——— = 10 log R= =-20 log R , (4-@) 
ref 


in dH re the distance reference, taken as 1 om. 


The spreading scheme for 4 refractive medium is shawn 
wei euire 4—-SC2]. This 1s based on ray theory which 
assumes that acoustic eneraqy does not cross the rays, with 
energy contained between two rays being conserved. The 
intensity at the reference ranaqe between the two rays shown 


is: 


= tt ty—o)) 





At a horizontal distance ¢r meters from the source, the 


imtensity 1S: 








Figure 4-3 Spreading loss in refraction. (From Clay and 
Medwinl2]) 








FE = 
—— a — = ee (4-10) 
# Soir 27,r Arjsin 8, maya /ecos 8, 
The reiation between intensities becomes: 
aA4T c = 
eae hes “6s INS Ly ince [eros Gz . 
= - = : ‘ (4-11) 
277 [\r|sin 8, | r Ar [sin 8, | 
The loss due to spreading is: 
I B8cos @ 
ir Oo 
= A : (At ae) 
er r {sin 8, | 
In terms of pressure, I, = Pr“) oes he. = pee Pata: 
and therefore 
Bet P1-1 @ cos 4 
ee Se A | = (4-17) 
eS Giaic Plein | 
Do Pa G 1 
Since Gy Sneii’s lew, cos SS eis ee yah 


~ 
to 


ae ; —P ne ees Oy 7 a pate 
pene a= r (\r|sin 8, | r Ar [tan 8, | : 


since en Oe —ermMesceaiaten bo aim 
excellent approximatian. 


Applying the dipole modei to this spreading loss 


eauatian gives 
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Equation 4-15 assumes a unique refractive path between 
source and hydrophone. In the Arctic there may be other 
paths due to the non-specular scattering of rays off the 
1ce canopy. Figure 4-4 illustrates how rays normally 
trapped in a surface duct, may be deflected down to a 
hydrophone. The minimum vertexing angle caiculated from 
the linearized sound velocity profile of Chapter 32 was 
m.0454 radians. For a ray to stay in a Surface duct above 
tre hydrophone it must be reflected from a slope of less 
BpoapeOres2 radians or about 2°. Tt 1s reasonable to assume 
that the ice canopy lacks local levelness to this order, s0 
that mon-specular rays must be accounted for. 

The rays which rebound from the ice canopy experience 
some loss. The attenuation for the FRAM IV experiment Ras 


been reported at 0.1 dB/km at 80 Hzlild. This attentuation 


may be converted to a loss per bounce. 


D 
GB = O.1 dB/km=_ , Ca= fey 
K 
where b = loss per bounce, and 
X = cycle distance. 


The cycle distance depends on the launch angle and the 
sound speed gradient. For a launch angle of O.O32 radians 
and the assumed sound velocity profiiese of Chapter =. the 


cycle distance 1s 3.7 km. Therefore, the loss per bounce 


is about 9.4 dk. This loss is low enmouach that even a ray 
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Figure 4-4 Specular and non-specular reflections frem the 
lce/water interface. 





which has bounced several times may contribute a 
Significant amount of eneray at the hydrophone. The non- 
specular rays cannot be ignored. Hounce loss at 
frequencies iess than 80 Hz are even smaller, since the 
data show a roughly linear dependence on frequency. 

To account for the non-specular rays the spreading 
loss is calculated using the ray averaging techniquelsSil. 
TRe oressure from 4a particular ray at a Given depth ana 


horizontal range, assuming a dipole source model, is 


AS sin”6, d&8 dr 


“(r 2) = cee 
a ; rftan a, | de X/2 


(4=37) 


dre 
The term ae represents the probability that a ray Bundie 
Will cross a certain depth, as shown in Figttre 4-35. For 4 


Single linear sound speed qradient the cycie distance can 


be written as 





> —_- : os eet 5 
sin = <r Sin ae aaa o (4-18: 


eee eo eee aS oF 
curvature, to a good approximation constant 
for all small angle rays ina linear sound 


speed gradient. 


Applying equation 4-18 to equation 4-17. and usina the 


small angl®= approximation, gives 


A= [ao] de, 
ieee = a (4-419) 


ie SB, 








Figure 4-5 Frobability of a ray bundle crossing a certain 
depth at a given horizontal range. 





In order to average the contributions of the possible 
rays, this pressure is intearated over all possible angles 
for a given receiver depth, and then averaged over depth 


down to the hydrophone at z 





o° 
ZG 2 
sy 2 A= 1 i 
Bo (is) = 22 dz J : (4-20) 
Oe ee (8, | 
O 8 


m 


where @, is the maximum launch 
angle of aray that wiil hit the hydrophane 
at a given range, and @, is the minimum 
launch angle. 


The angle 8, is a function of 8, and z 
e 


z=er [cos 8, - cos ees : (4-21) 


Using the small angle approximation for cosine leads to 





a 
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Substituting this into equation 4-20 and evaluating the 


integral over angle gives 





aan 
“ — A =“ 
07 (a dz eta — 2S) 
r~ re “65 
> <2 
6, = — for all = so the second term within the integral 
iz 
Cc 
1S always zero. Evaluating the first term aver centh cives 


Seer eesoiom for pressure in terms at r and @,. 
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a A= a et 
po~(r)} = Y-" -—- {[8&*- - ns (4-24) 
Si y ¥ i" 
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This expression can be used to find the source strength. 


: (4-23) 





The spreading function, G, 15 Gresented as a function Gf 
alone since Sy depends on r. For each r there is a unigue 


8 and therefore, a unique G. The spreading function was 


wi 
calculated for horizontal ranges from 300 m to 20,000 om, 
and tablulated in Appendix D. -The spreading function is 
shown in a log-log plot in Figure 4-6. For comparison the 
equivalent spherical spreading for a dipole source 2 alsa 
shown. From this one can see that source strenaths 
Calculated using the spherical soreadina taw iead ta an 
unrealistic cependence con rande. 


Mpeetes FO get © and r- a iinear sound speed 


Cc 
gradient of 9.954 sec! was chosen. This gradient gives 
eaves eae © — S&S km as the multiple step profile used 


in Chapter 3. Three kilometers was chosen since it was the 
median horizontal range for the noise events. 
The spreading loess function and measured seak oressure 


maanitudes were used to caiculate Ginols strencth. 


ie 


F, = 2AA = 2ApoV/Gir) - aa 
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Figure 4-6 Spherical (heavy line) and refractive spreading 
loss, 


Be), a5 a function of horizontal range. 





For a particular event the scurce strenath was 
calculated from each hydrosnone, and then the average taken 
as the source strength for the event. The standard 
deviati an within each event ran from 9% to 125% of the mean 
Value. An event strength summary in Appendix C lists the 
mean measured peak pressure and the mean dipole strenath, 
along with the standard deviation on these values for each 


event. 


strength of Background Noise 

IT was interested in the effect that environmental 
loading might Rave had on the temporal, spatial and 
strength statistics. It has been shown by Makris and 
DyerlC1i0O] that low frequency (10-20 H2 band) ambient noice 
rms pressure, averaged over a long time, correlates weli 
with environmental stresses and moments. Since [ haa 
ambient noise pressure for most of the period of the FRAN 
IV experiment, and since [ Rad environmental stresses and 
moments available for only a part of the time, [ chose to 
use the 20-80 Hz lona-time-average rms pressure as my 
environmental indicator. 

The 10 to 20 Hz band ambient noise pressure was 
converted to 20 toe 389 Hz band pressure in the folilowine 


= 


Manner. Figure 4-7 shows the typitcai spectrum for centrai 


Arctic pack ice noise. The portion of the spectrum between 


19 and 100 Hz can be approximated By a straight line. 
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Figure 4-7 Composite central Arctic ambient noise spectrum 
observed during the FRAM IV experiment. ‘From Dvyerlé]) 





Fog 7s =A — lag tJ] + BE , (4-27) 
where A = sioace = -1.7275 FPa*/Hz* = 
B= intercept = -1.0909 Pa~/Hz ‘ 
or 
S=10F.4 _ (4-28) 


Tre band rms pressure relates to the esnoectrai level oy: 


z - =e 
Beenees. = Js df . (4-29) 
b 


T have assumed that as the sound pressure level changes 
from time to time the intercept ‘BH changes. but the slone 
remains the same. By substituting equation 4-28 into 4-29. 
and using the known 10 to 320 Hz ambient noise band, B can 


Be weritten in terms of the known oressure. 


oa 


P'rms. 10-20 


SS Sos (4-39) 
Ky 
where 
De eat 1) 
Ky = —_—_—__—______—_| = 0.1020 Hz . 
A+ 1 


The ambient noise rms oressure for the 2060 to 80 Hz Gand may 


now be found. 


BO 
= ay ~HeA —=41} 
Pp rms.29-B80 — f Lee df . (4— 2.1) 
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hen 
a 2 
=] 220s. lo—26 Ky : 


where 

; go tA ot a ee 29 6A tie) | 

a $$ queii__ d = «SCF 0988 Hz ; 
= A + 1 


Gm finally. 


: (4-32 


Prms.20-80 ~ Prms,10—-20 





Prms,2z0-80 — 9-78 Perms ,io-zo ° 


Thus the band from 20 to 80 Hz 16 virtually identical to 
Siesvaesrrom 10 to <0 He in rms pressure, for long—-time— 
averages, and in turn, is an acceptable surrogate for 
environmental forcing (annlied stresses and moments). TRe 
=9 to 8@ Hz Band ambient noise rms pressure for each at the 


tapes investigated is found in Table 4-1. 





Tabie 4-1 =O to BO Hz Hand Ambient Noise rms Fressure 


Tape # Date Recorded Perms, 20-80 
(Fail 
4001 na ge eee Not Avaiiabie 
~OO1 ey oy ce Not Available 
2009 3-30-82 ©. O22 
S001 5-31-82 QO.019 
4003 4-01-82 0.044 
4005 Zot oe 1B) UE a 
4007 Seo Ce Oe. 
4009 fe 0.010 
4011 4-02-82 0o.010 
4013 4-O2-82 Oa 135 
2023 4-08-82 0.037 
4015 4-09-87 O. O40 
5047 At S5—o eS 0.010 
4016 Bao ee O.017 
4019 4—S—ae O.016 
4021 4-19-82 O,Q13 
4023 S19 ee O.O1l 
4024 4-19-82 er oC) tal 
4027 ee O.0O12 
4029 4—7O-—82 MeAg tS 
4031 §—29-82 a 
4033 4—7?0-82 0.012 
4046 4—255—-38 > G,. 034 
4047 Fela O.114 
4049 A= 2A © 0.140 
4051 4 O.140 
40355 ey Pe 0ao 
4055 Ae 0.982 
4057 B= 72-82 OS 
4059 q—2a2-82 0.065 
4061 2 ae 0.034 
4063 ceo oe 0.0293 
4oG65 oe (yee 
41567 qe eS ee 





CHAPTER 5S 
ANALYSIS QF NOISE EVENTS 
Detection Analysis 

A total of 34: tapes was examined, for a total time of 
662 minutes. (For a few of these tapes the entire 20 
minutes was not used.) 

There was a total of 499 detections of events flagged 
om at least SO% of the hydrophone channels. Of these, 139 
were man-made artifacts, and 125 were false alarms 
(detections which were so weak that no pattern for taking 
time delays could be discerned). There were 199 unique 
events, and 36 multiple dtections of those events. Stated 
in another way, of the detections which were not artifacts, 
63.3% were strong enough to support analysis and 34.72% were 
too weak to reasonably analyze, and hence labeled false 
alarms. 

Since the detection process depends on signal-to-noise 
ratio, the level of background ambient noise should affect 
the event detection rate. Figure 3-1 shows normalized 
ambient noise pressure, number of false alarms per tape, 
and number of unique events per tape for each tape 
examined. There is some trend for more events being found 
when the ambrent pressure is lon, and more false alarms 
declared when the ambient pressure is high. 

This is more clearly seen in Figure 3-2, which shows 


the average number of false alarms and unique events found 
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S-i Normalized ambient noise pressure, number of 
events per tape, and number of false alarms per 


tape for each data tape examined. 





pressure. 


Events = \\\) False Alarms 





per tape in each of four background noise pressure ranges. 
The 0.01-9.02 Fa range used 15 tapes to compute its 
average, the 0.02-0.03 Pa range 4 tapes, the 0.03-0.04 Pa 
range S tapes, and the over 0.04 Pa range 8 tapes. Two 
tapes were recorded during the first few days of the FRAM 
IV experiemnt, before the 10-20 Hz band ambient noise 
recordings were started. 


A Breakdown of detections for each tape is found in 


Appendix C. 


Temporal Analysis 


The interarrival time between events ranged from 1 to 
1064 seconds. Each event time was taken to the nearest 
second, and no events were taken as Raving the same event 
time. If two events happened in the same second, ane was 
judged to be earlier, and the two events were given event 
times one second apart. The interarrival time for a 
particular event was measured from the previous event, 
except for the first event of a tape, which was measured 
from the start of the tape. 

The interarrival times were divided into bins of 20 
seconds. The first bin ("0") contained events which had 
interarrival times from O to 19 seconds, the second bin 
from 29 to 29 seconds, and so on. The number of events per 
bin is presented in Table S-1 and shown graphically in 


Figure 5-3. A complete listing of interarrival times for 
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Nurnber of Events 





20 Second Bin (0=0-19s, 1=20-39s, etc) 


Figure S-S Number of events found per interarrival time 
bin. 





each event is found in Appendix C. The mean of the 
interarrival times CLL ) 1s 100 seconds, and the standard 
deviation (Q) 166 seconds. In terms of bins, the mean is 


3 and the standard deviation 8. The standard deviation 15 


1.66 times the mean. 


Table S-1 Number of Events per Interarrival Time Bin 


Bin Events Bin Events 
O 6&7 
1 36 1 O 
Re 19 Ks O 
3S 15 33 1 
~ 1Z 34 O 
2 7 SS O 
S 4 36 O 
Z 3 37 O 
8 S 38 O 
9 1 39 0 
10 S 40 O 
ii 4 4i 1 
12 2 42 0 
13 1 43 O 
14 3S 4 0 
15 Z 45 O 
14 0 46 O 
17 0 47 O 
18 O 48 0 
LS O 49 1 
20 1 50 O 
an 0 31 0 
ZZ 0 a2 0 
Za Z 93 1 
24 1 54 0 
Za 1 2 O 
26 1 5&4 0 
27 Z 7 O 
28 1 98 0 
29 1 59 O 
3O O 





Three different probability density functions were 
investigated to find an appropriate fit for Figure 3-3. 
They were 1) a half-gaussian distribution, 2) an 
exponential distribution and 3) a J shaped distribution. 


The half-gaussian probability density function isf3i1: 


2 ? 
-_ cen ee a 


o(t) = —— e : (S-1) 


The general equations for mean, mean square value and 
variance (O67) Can be used to solve for the unkown 


constant, t,: 


CO 
jo -/ BiG ee a (Sq) 
QO 
CO 
mean square value -[ mCeE) Sac. . (a=) 
OQ 
eco 
g 2 =f ce - Mo? pee dt. (5-4) 
OQ 
= mean square value — ju * : 


Substituting equation S-1 into equations 3-2, 5-3 and 5-4 


leads to the following relations: 


Oo ? , Oo 


27 
t. = Le ; mean square value = t = Ca a) 


@ + = “ae to- 3: OF 0.756 { - 





This value for t, was used in equation S-1, and the 
probability density function integrated over appropriate . 
limits to get the number of events in each 20 second bin. 
The result 156 Sat ced against the experimental distribution 
in Figure 3-4. 

The second distribution (the exponential) belongs to 
the family of gamma distribution functions [15]: 


1 OQ -t/t, 
ee tl : (5-6) 


Pee lacy 


When (Y = 0, this becomes the exponential probability 


density function 
: CS—7) 
= — Fe cm 
P © 


Cc.” 


Again using equations 5-2, 3-3 and 3-4 leads to: 


eh =t, :; mean square value = eile : (5-8) 


oe O= pf . 


The exponential probability density function was integrated 
Over the bins, and the results are shown in Figure 5-5. 
Another demonstration of the fit of the exponential 
probabilty distribution is shown in Figure 5-64. Taking the 
natural log of the function should lead to a straight line 


when plotted against time or bin number. The straight line 
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Figure 3-4 Half-gaussian distribution compared with 
experimental values. 





70 
60 
50 
2 
c 
; 
4 4) 
0 
Ws 
mh 30 
C 
5 
Zz 
20 
10 
0 = aS ees 


0 2 4 6 8 10 12 +14 «+16 ~~ 18 20 
20 Second Bin (0=0-19s, 1=20-J9s, etc) 


Figure S-a Exponential distribution compared with 
experimental values. 
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Figure S-& Semi log plot of exponential distribution and 
experimental values against bin number. 





in Figure 5-6 is a plot of the natural lag of the points 
calculated using the exponential probability density 
function. The experimental points seem to curve rather 
than 11e€ on a straight line. 

The mee distribution (J shaped) is also a gamma 
distribution. The J shaped distributions are characterized 


by CY < 0. $I chose a fairly common distribution with 





CY = -9.5 . The probability density function is: 
i lee ta ae 
pe t e ; (5-9) 
V toll 


and the key parameters are: 





= =U : mean square value = a ; CS fo) 


This distribution is plotted against the experimental 
values in Figure 3-7. The natural log of both calculated 
and experimental points are plotted against bin number in 
Figure 3-8. This distribution seems to fit the 
experimental points best of all. The J shaped probabilty 
density function goes to infinity at zero, but 1t is 
imnzearable. 

A Chi square goodness of fit test was done on all 


tnree distributions. The results are summarized in Table 
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Figure S-7 J shaped distribution compared with 
experimental values. 
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Figure 3-G Semi log plot of J Shaped distribution and 
experimental values against bin number. 





Doe Also presented in Table 5-2 are the ratios of 


standard deviation to mean. 


Table S-S Comparison of Distribution Functions 


Chi square OU 
Experimental = ,20oG 
Half—-gaussi an P27 ee 0.76 
Exponential 45.69 ioe 
J shaped 10.65 1.41 


For a distribution to pass a goodness of fit test it 
must have a Chi square less than a prescribed limit. The 
limit for my test (9 degress of freedom, CY = 0.005) was 
25.6016]. Only the J shaped distribution passed the Chi 
SQuare test. It also has CAG. closest to the 
experimental values. In summary, the interarrival data 
reasonably fit a J shaped distribution given by: 


; -1/2 -t/2 
pt) = ot e : (5-11) 


4/2 TTL 
Since event detection rate depended on ambient noise 
level, interarrival time between events should also show 
environmental dependence. Table 3.3 gives average and 
standard deviation of the interarrival time for different 


ambient noise pressure levels. 





Table S-3 Background Noise Level Dependence of 
Interarrival Time 


Ambient Noise Interarrival Time 
rms pressure 
(20-80 Hz) mean standard deviation 

(Pa) (sec) (sec) 

O.01-0.02 &7 129 

0.02-0.03 190 153 

Q0.02-0.04 147 143 

over 0.04 183 318 


The tapes having a background noise level of 0.01 to 90.02 
Pa have a significantly shorter interarrival time than 
tapes in the other three pressure groups. As with 
detection rate, the interarrival time does depend on 


ambient noise level. 


epatial Analysis 

After removing nonlocatable events and events located 
outside a horizontal range of 20,000 meters, 164 events 
remained. These were grouped by horizontal range into 42 
annuli of equal area as shown in Figure 3-9. Each annulus 
1s a 350 Square km ring centered at the array origin. The 
first annulus ("0") went from O to 3090 meters, the second 
from 2090 to 4270 meters, and so on. 

Table 3-4 shows the number of events per annulus and 


Figure S-19 shows this distribution graphically. 
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Figure 5-10 Number of events per radius annulus. 





Table 3-4 Number of Events per Annulus 


Annulus Events Annulus Evants 
O ot zl O) 
1 We “2 0 
= + ae il 
os a =4 i) 
+ ee el a 
= = =O 1 
5 9 ai QO 
eo a eS O 
8 1 ay, CG 
9 = td O 

10 = 1 aj 
11 a ee O 
ee + ao 1 
a 1 34 O 
14 1 Bs 1 
1s + ma i) 
146 Q = 1 
17 oO E 1 
18 1 & O 
19 © 40 ij 
20 Cc 41 1 


The average number of events per annulus is 3.93 and the 
Standard deviation is 14.135 events. Figure 3-190 shows that 
the number oF events Tound is highly dependent on their 
range from the array. In the center annulus there were 
over 20 times the mean number of events. 

The dependence on range is not a surprise, since 
spreading (and possibly scattering and other losses) wili 
reduce the strength of weak transients down to the ambient 


moise level. For this reason, the center annulus is 


> 
ts 


R= 


sircbDably tne cest indicator oF actual event density. 


tnis ring there were 91 events per 20 square kilometers per 


‘| 


662 minutes of observation or aporeximately @.4 events per 





squaré Kifiometer oer hour. 

The averace number of events per annulus and the 
nunber of events per square kilometer per hour should 
Cepend on Background noise level. The average number oat 
events per annulus was found for each ambient noise rms 
pressure range, and adjusted to reflect the number of 
events in a 4662 minute period. The results are seen in 
Tasle 3-3. The number of events per square kilometer ser 
nour for the center annulus are also shown in Table 3-2. 


Table 3S-S Average Number of Events per Annulus 
for 4 Ambient Noise Levels 


Ambient Noise Events Minutes Adjusted ae te 
rms pressure per of tape Events oer per km* 
beo—GO Hz} Annulus Examined Annulus per Ar 
Oeel—O.02 Fa =. 60 2 ac, a Sieh tie 0.452 
Seee-O,05 Pa Oaee 7 7a 2.46 0.250 
GeO. 084 ra 0.45 i199 emi? Geto 
over 0.04 Fa GaaiG dep 1230 2076 

Emherre 
Fopulation a. 7S 662 ee yas Onela 


The average number of events per annulus and the 
maumber of events per square Kilometer per hour beth retlect 


the ettect of signal—-to-norse ratio an the detection 


The entire population af Svents was investigated far 


angular denrendence. Figures 2-11 and 3-12 show the number 
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Ficure 5-12 Number of events per 70° sector. Radius 
Gives the number of events, while angle indicates the 
sector measured from the northern leg of the array. Each 
ring represents i909 events. 





of events found per 20° sector. In the polar diagram 
(Figure S-iS) the radius shows the number of events. The 
angles are measured from the northern leg of the array. 
Figure S-13 is a polar plot showing the number of events 
per 10° sector. There was no predominant angular dzrrectian 
found. However, some preference can be seen for bearings 


of 320° and 190° from the northern leq of the array. 


Strength Analysis 

The mean hydrophone peak pressure magnitude for each 
event feil within a fairly narrow band of values. The mean 
peak pressures ranged from 1.32 to 0.16 Fa, with an average 
of 9.356 Fa and a standard deviation of 0.20 Fa. Figure 3-14 
shows the mean hydrophone peak pressure values for ail 
events located between 190 m and 20,000 m plotted against 
mange from the array origin. 

The different symbols shown in Figure 3-14 represent 
events dGuring each of the four ambient pressure categories. 
The events with a Arq@her mean hydroohone peak oressure have 
a tendency to accur during higher ambrent norse levels. 
This can be seen in Table 3-6, where the maximum, minimum, 


average amd standard deviation of the mean hydrophone peak 


in 


pressure values are given for each of the four ambient 
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Figure S-isS Number of events per 10° sector. Radius 
gives the number of events, while angle indicates the 
sector measured from the northern leg of the array. 
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Figure 3-14 Mean hydrophone peak pressure measured for 
events Between 100m and 20,000 m, plotted against 
horizontal range from the FRAM IV array origin. 
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Table 3-6 Mean Hydrophone Feak Fressure 
for 4 Ambient Noise Levels 


Ambient Noise Mean Hydrophone Feak Pressure 
rms pressure (Pa) 

WaeeeQ Hz) MAN min average std dev 
Pemoi-O.02 Fa Lee O.16 eer ie We 
Peee—-O.0S Pa ve ae Of 22 Ono eee) 
Gee. 04 Fa eve ey al eas O. 43 er ad 
over 0.04 Fa ie cs O. 16 O.49 OG. 24 
Entire 

Fopulation =e 0.16 0.36 QO. 20 


Source strength (Fi? was found for the events which 
had Hydrophone locations between 3200 m and 20,000 m from 
the event. 

The dipole strengths ranged from ZZ kN to 4.9 MN, 
With an average of 431 EN and a standard deviation of S332 
mM. The distribution of strengths for the isi events 
evaluated 1s shown in Figure 3-135 and in fable s-?7. 

Figure 3-146 shows the dipole strength for all events 
plotted against horizontal range from the array origin. 
Again, tt can be seen that the stronger events accur when 
the ambient pressure level is high. Table 3-8 gives the 


Sstrenath values for the different ambient noise leveis. 
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Distribution for a Population of Events 


KN 
kN 
kN 
kN 
kN 
kN 
kN 
kN 
KN 
kN 
kN 
kN 
kN 
kN 
kN 
MN 
MN 


# of Events 


7, 
Za 
33 
23 


NerFP REANIM OP on 


Strength versus Ambient Noise Levels 
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Figure S-16 Dipole strength for events between 200 m and 
20,000 m, plotted against horizontal range from the array 


Ome Qin. 
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Most of the events evaluated for strenath occurred 
during the lowest ambient noise levels. Figure 5-17 shows 
the strength of events that occurred when the ambient noise 
was 90.01 to 0.02 Pa. The log of the dipole strength is 
plotted against ane log of the horizontal range from the 
center of the array. The points scatter more so to the 
upper left rather than lower right, because distance itself 
fiiters out weak events. A weak signal from far away would 
not reach the Rydrophone array with enough amplitude to be 
distinguished from the background noise. And events 
located farther away would tend to be strong events. 
However, events located close to the array should have the 
entire range of source strength levels. This would proecuce 
a wedged shaped plot of weaker events close to the array. 
Indeed, Figure 3-17 shows a general scattering with perhaps 
a wedge of weaker events near the array oridin. 

Nonetheless the trend shown in Figure 3-17 suagests 
that the ray average mogel used to estimate refractive- 
surface reflective spreading may need to be replaced with a 
more refined model. For example, horizontal ranaces iess 
than about 1000 m mav include too small a loss, and 
therefore lead to too small a strength, because the 
reflective contributions may not be as large as imputed. 
Sucn & criticism is supoeorted by the notion that for a 
Geen steppe. € 5 Pertective rays and hence ray averaging 


ccecurs only beyond a critical horizontal range. 
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Figure 5-17 10 log of Dipole strength (dB re 1N) versus 
19 log of horizontal range from the array origin (dB re 1 m) 


for events occurring during an ambient noise level of 0.91 
ee O.02 Fa. 
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The foregoing speculation suggests that the average 
dipole strength for the lowest ambient noise case is best 


found from the events farther from the array, and is 


F, = 10° N™ 320 kN, ae 


with a much smaller standard deviation than in Table 3-8. 
Presumably corresponding adjustments could be made for the 
higher ambient noise cases, but the FRAM IV data set 
contains too few events at higher ambient noise to plot as 
Pee gure 3-17. 

The strength analysis 1S a somewhat ambivalent one 
because of spreading model uncertainty, and because data on 
ice slopes are not available. But the dipole picture of an. 
event likely has some validity, and at least rough 
estimates of its strength Rave been extracted from me 


data. 
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CHAPTER & 


SUMMARY AND THOUGHTS 


Through the use of a detection program, visual 
confirmation and a location program, a population at 177 
Arctic noise transients was gathered. There are four major 


Besuits. 


ie 


First. mare events are found when the ambisnt oressure 
is low, ard more faise alarms when the ambient pressures is 
hich. The interarrival time and the average number of 
events per unit area aiso depend on ambient noise level. 
Since more events are found when the ambient noise is low. 
the interarrival time decreases, and the spatiai density 


increases. 


Second, the interarrival times were fit to several 


ui 


possible orobability distributians. The interarrival time 
distribution best fits a J shaped gamma distributison. phe 
mean imterarrival time is 1990 seconds. 

Third, the number of events per unit area 15 nighly 
dependent an range, since distance filters out weak 
transients. The event density in the annulus closest tc 


— 
— 


the center of the array was 9.53 events per square kilometer 


Nt 


cer Four over ali observations and 9.3 events per sauare 


tl 


OME : 


pend 


hi r Ger hour ror quiet times. There 15 no 


! 


predominant ancular dependence to the spatial distributian 


9 events. 
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Last, the mean dipoie strenath for the observed events 
25 4430 KEN overail and 249 EN during Low ambient noise 
ieveis. stronger events occurred during high ambient noise 
levels. A refinement of the spreading loss model used to 
Calcuiate these values may lead to values whicn are 


Slightly higher. 
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is a GBreric acoustic Sevents is far from 


ap 
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a 
Ud 
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thi 


complete. several areas for improvement have been 
mentioned earlier in the thesis. The detection program 
needs to be made more robust to eliminate the event time 
error. A scheme for ignoring artifacts should be includec. 
The loecation orogram wastes time looking in the wrong 
direction, aithough the Bearing accuracy of the program is 
very qood. The aigorithm should be changed to quickly fine 
the right Bearing, and then search in a sector. 

The type of Beach event, whether it was a pop aor 4 


imtormatiaqn and 


C. 


Whine, eas not recorded. Peete cerme thi 


rm 
ift 


correlating it with interarrival time and range still needs 


to be done. 
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APPENDIX A 


User’s Guide for the Adetect Frogram 


Figure A-1t: Flow Chart of the hdetect Frogram 


source Code for the Adetect Frogram 





See io 


USER’S GUIDE FOR THE Adetect PROGRAM 


The purpose of the Adetect program is to detect 
ambient noise transients amidst the background ambient 


noise recorded on a FRAM data tape. This is done by 


t 


comparing the short average of data points to the long 


average of points on a single channel in order to flaqg a 


Bossibie cgetection, and then waiting until 3SO% of the 


iF 


crannels are flagged to declare an actual detection. 

The input for the Adetect program is a framread output 
file without headers. A FRAM data tape is read into the 
file by the command 

-framread —-head neni 


The program will ask for the input device (tape drive 


designation), the output file, and the number af data 
seqments to skip and to read. Fach seaqment represents =.6 
seconds of data on 24 channels. The tremread program reads 


a first seament which contains no data records, s0 you 
should specify skipping one more segment than you would 
normally calculate. For example, reading the entire first 
half of a 20 minute FRAM IV tape would require the response 
of 
1 160 RETURN > 

£G tre aquesticn of “enter #ekip, #segqments:". 

Once this input file has Been created the Adetect 


orsearam can be used. The program is started with the 





alo 


command 

hdetect RETURN > 
The program will ask for the FRAM tape number, the Julian 
date of the tape, and the start time of the tape in hours, 
minutes and seconds. The program will then ask you to 
select the channels you wish to use. In most cases the 
FRAM data tapes did not have ambient noise hydrophones tied 
into ali channeis, and the specific channel that a 
hydrophone was recorded on changed throughout the 
experiment. Which channels were in use and for which 
hydrophones can be found in the experiment logs. The 
program assumes that the channel number 15 equal to the 
hydrophone number, but allows you. to change this by 
inputting the channel number and the proper hydrophone 
number, ar "O" if the channel is not in use. For example, 
if channel 3 was not used , and channel 7 was used for 
hydrophone 21, the input would be 
“RETURN > 


Be <RETURN> 
, <RETURN> 


On 
On a2 


The "0,90" ends the changes to the channel selection. You 
must now hit any key to continue the program. 

You will be asked to enter the input device (the input 
framread file), the number of skips and segments, and the 
neme of the output file. The output file does not have to 
exist before the program 1s started. It will be created by 


the program. The number of skips and segments are those 





Sa 


that you would calculate using 3.8 seconds per segment. 
For example, to proces the entire first half of a 20 minute 
FRAM IV tape the number of skips and segments would be 
Q,160 “RETURN > 
This is all of the input required by the user. The 
program proceeds from this point without user interaction. 
The output of the Adetect program is a file containing 
a list of detections in the following format: 
tapenumber Juliandate hour minute seconds 
O eventnumber eventtime 
Channel hydrophone timedelay amplitude 


channel hydrophone timedelay amplitude 
Channel hydrophone timedelay amplitude 


channel hydrophone timedelay amplitude 
QO eventnumber eventtime 
Channel hydrophone timedelay amplitude 


channel hydrophone timedelay amplitude 
QO eventnumber eventtime 


mp 
The "O" at the start of a line indicates a new event 
detection, and the "-1" at the start of a line indicates an 
end of file. Each channel that was flagged for a 
particular event is listed with its hydrophone number, 
timedelay from the earliest channel signal arrival, and its 
peak voltage amplitude. This outfile can be used as the 


imput file for the location programs without modification. 
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/* startdoc 
hdetect.c 


program to read whoi-segy format data tapes from the fram IV program. 


usage: 
hndetect 


program is interactive. 


by Mary Townsend-Manning 


enddoc 
a] 


#include <stdio.h> 

#include <math.h> 

#define NCHAN 25 

#define RECLN 950 /* number of samples per trace (4 bytes per sample) */ 
#define OBYTES 3800 /* number of bytes per record output */ 
#define ZERO 0 

#define LONGFILTLN 64 /* length of long average filter */ 
#define FLN 64 

#define MAXEVENTS 4 

#define RESET_DELAY 0.3 

#define RATIO 2.38 

#define SIG DELAY 0.02 

#define THRESHOLD 0.5 

#define EVENT_DELAY 0.5 

#define END -1l 


double vconv(x,y,n) 
vregiceer float *x, *V-; 
EeQUSeer Int <n; 


register double sum = 0.; 


rf 2 0) 
{ 
eo. { 
Sie et as 
\ewnaley==ne > 0); 
} 


return (sum); 


main () 


Source cade for the hdetect program. 
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float l_ave(NCHAN], sh_ave, timeofflag(NCHAN]; 

float ampofflag(NCHAN], firsttime, timedelay (NCHAN] ; 
float chaneventtime (MAXEVENTS] (NCHAN]; 

float chaneventamp (MAXEVENTS] (NCHAN], eventime ([MAXEVENTS] ; 
int flag (NCHAN]; 

int num_active_events; 

anc em, et: 

int event_flag(MAXEVENTS] (NCHAN], flag sum, number _of_events; 
int event_number (MAXEVENTS] ; 

int nchan = 0; 

iitwchariuNGHAN|, kK, 3, 2, tid, date, hour, min, sec, 1; 
Ent channel, data, toggle; 

char answer; 

float A[(FLN), longfilt (LONGFILTLN] ; 

int nskip,nseg ; 

char oddobuf (NCHAN] (OBYTES],evenobuf (NCHAN]) (OBYTES]; 
int €ount, error; 

char fname(80)],iname(80]; 

FILE Bt c, Pee, ~tp, *fopen () ; 

float timeseries (NCHAN] [RECLN+2*FLN]; 

float time = 0.; 


num_active_ events = 0; 


for (i1=1;1<25;i++) { 
enan(ij=(2) ; 
flag (4)* 0; 
for (m=1;m<MAXEVENTS ; m++) { 
event_flag(m] [iJ = 0; 
chaneventtime([(m] [iJ = 0; 
chaneventamp(m] (i) = 0; 
} 
} 


/* Program initialization from keyboard */ 


Eprine. (stderr, Program Initialization\n™) ; 
fprintft(stderr,"enter FRAM tape #\n"); 
EPscanttsctdin, td’ ,<tid)-> 


fprintf (stderr,"enter Julian date\n”); 
f£scanf (stdin, “td", édate) ; 

fprintft (stderr,”"enter time - HR,MN,SC\n"); 
EScane (Stain, ta, $a, $a" ,anour, &min, &sec) ; 


fprintf (stderr,"default values for channels and phones\n"); 
fprintf (stderr,”"are channel # = phone #.\n"); 


fprintf (stderr,"enter channel, phone to change.\n"); 
fprintf(stderr, ”"enter ‘0’ for phone, to eliminate a channel.\n"); 
Poei nce (stderr, enters’ 0,0 to quit. \n"); 

fscant (stdin, "td, $d", é), &k); 


while (3 != 0 &6& j < 25) { 
chan( 3) = 7k; 


Scurce code for the hdetect program. 





£Scant (stdin, %+d, sd°,&}, &k) : 
} 


fprintf(stderr, "FRAM TAPE %d Julvan Date: td\n", tid, date); 


EPEIntE (stderr, “ Time: tata 2td\n", hour,min, sec); 
for(iwl] -a<l3eats) 
ferantzistcerr, "CH Bd PH <td CH %d PH $&d\n", 


2,CNan li) 7itiz,enan(it1Zi)- 
/* Check to make sure inputs are correct -- Change if necessary */ 


scant (stdin, $c”, answer); 
fprintf(stderr, "Hit any key and RETURN, when ready."™); 
fscanf (stdin, “tc", Sanswer); 


for(iz1l;i<NCHAN;i++) { 
afi(chan(ij) != 0) nchan++; 
} 


forint’ (stderr awenter input device: “); 
scanf ("%s",iname) ; 


if((iptr = fopen(iname, "“r")) == NULL) 
{ 
formtt (swderr,"can’t open %s\n", iname) ; 
exit (1); 
} 


forintf (stderr,"enter #skip, #segqments: \n"); 
fprintf(stderr, "values of 0 and 320 will read entire tape\n"); 
scanf ("%d, td", &nskip, &nseg) ; 


/* load bandwidth filter */ 


if((fp=fopen("“"PMf£loat","*r")) == NULL) { 
printf ("cannot open bandwidth filter file\n"); 
exat (0): 
} 
for (i=0;1<64;1i++) 
Eacane (ep, si ,et (i) ) > 
felose (fp); 


/* load averaging filter */ 


for (1™0; i1<LONGFILTLN; i++) 
longfilt{({ijJ = 1.0/ (float) LONGFILTLN; 


/* Open output file */ 


Source code for the hdetect program. 
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foErnee (secerz, enter out=file; ™"); 

scant ("%s", fname) ; 

if( (ptr = fopen(fname, "w")) == NULL) 
{ : 
fprintf(stderr, "can’t open %s\n", fname}; 
exe (1); 
} 


fprant£ (ptr, "td td td td) td\n", tid, date, hour, min, sec); 
time = time + 3.8 * (nskip-l) ; 
/* enter first record */ 


if (nskip%2 == 1) { 
for (j=l; j<NCHAN; j++) ( 
fread (&evenobuf[(j](0], sizeof(float), RECLN, iptr); 
} 
toggle fl; 
} 
else { 
for(j=1;j53<NCHAN; j++) { 
fread (aoddobuf[j](0], sizeof(float), RECLN, iptr); 
} 
toggle = 0; 
} 


/* ENTERING RECORD READING MODULE */ 
time = time + 0.504; 


/* forinet (scder=, Using butter size $d bytes\n", sizeof (buf)); */ 
for(rel?: 2 < nseg;: i++) 
{ 
7* fprineS (stderr, time = Sf \n", time)? 
GOcInce (stacrr, —precessing cecora td\n°, nskipti); */ 


/* read next record into appropriate buffer */ 


if (toggle == 1) { 
for(j=1; 3<NCHAN; 53++) [ 
fread (socdcbesi 7) (0), sizeor(tlicat), RECLN, iptr); 
} 
toggle = 0; 
} 


else { 
for (31; <NCHAN; j++) { 
fread (&evenobuf(j)(0], sizeof(float), RECLN, iptr); 
} 
toggle = 1; 
} 


7/*-filter and square data */ 


Source code for the Adetect proaram. 





Lf (toggle mari). { 
for (j=l; j<NCHAN; j++) { 
if(chan(j) != 0) { 


sq_filt (&oddobuf (j) (0), &evenobuf(j) (0],h, 
&timeseries(j] (0)); 


l_ave{j] = vconv(élongfilt[(0],&atimeseries(j] (2*(FLN-1)], 
LONGFILTLN) ; 


} 
} 


else { 


for (j=1; j<NCHAN; j++) { 
if(chan(j] != 0) { 


sq_filt (&evenobuf (j) (0), &oddobuf(j](0),h, 
&timeseries(j] (0]); 


l_ave(j]) = vconv(élongfilt (0), &timeseries(j] (2*(FLN-1)], 
LONGFILTLN) ; 


/* ENTERING EVENT DETECTION MODULE */ 
for (k#0;k<RECLN;k += 4) { 


for (l=1;1<NCHAN;1++) { 
Ff (chant 1) te" 0) 
l_ave({1] = (63.0*l_ave(1] + timeseries([1] (2* (FLN-1) +k]) /64.0; 
sh_ave = (timeseries (1) (2*(FLN-1)+k]) + 


timeseries (1) (2* (FLN-1)+k-1] + timeseries (1) (2* (FLN-1)+k-2] + 
timeseries (1) (2* (FLN-1) +k-3)) /4.0; 


J* zeset cold flags */ 


if (flag(1) == 1 && 


{time - timeofflag(1]) > RESET DELAY) 
flag(1] = 0; 


/* set flag if RATIO of signals is reached */ 


if ((sh_ave/l_ave(1])>=RATIO) { 
if (flag(1] == 1) { 


if (sh_ave > ampofflag(1l]) { 
timeofflag(1l]) = time; 
ampofflag(1] = sh_ave; 


} 
else { 


if (num_active_events == 0) { 
flag (lies i: 
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timeofflag({1] = time; 
ampofflag({1l] = sh_ave; 
} 
else { 
for vimtiem= (num active events, #1) sm+t) { 
LE C¢event_flag({m]) (1) == 1) { 
1£ {{time - chaneventtime([m)] [1)) 
<= SIG_DELAY) [{ 
if (sh_ave > chaneventamp[m] (1)) { 
chaneventtime(m](1] = time; 
chaneventamp(m] (1) = sh_ave; 
} 
} 
else { 
flag({1) = 1: 
timeofflag({1l] = time; 
ampofflag({1] = sh_ave; 
} 
} ‘ 
else { 


event_flag(m) (1) = 1; 
chaneventtime(m]({1]) = time; 
echaneventamp(m] (1) »“sh_ave-; 
m = num_active events; 


} 
/* end of set flag module */ 
/* start new event module */ 
flag_sum = 0; 


for (l=1;1<NCHAN;1++) { 
2£ (flag(1] "ssl) flag sumt+; 


if (((float) flag_sum/(float)nchan) >= THRESHOLD) { 
num active Tevenes +s. 
eventime [num_active_events] = time; 
number of events++; 
event_number (num_active_events] = number _of_events; 


for (l=1; l<NGHAN;1++) { 
event_flag(num_active_events] [1] = flag{1]; 
chaneventtime (num_active events] (1) = timeofflag(1); 
chaneventamp (num_active_ events] [1] = ampofflag[1l]; 
flag (Li, 0; 
timeofflag(1l] = 0; 
ampofflag(1) = 0; 
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} 
/* end of new event module */ 
/* start of deactivate old event module */ 
if (num_active_events > 0 && (time - eventime[1]) > EVENT_DELAY) 
forikt. (Otc, 6G sc. sc \n, » ZERO, 


event_number[1],eventime([1]); 
Vis fprintf (stderr,"td %f\n",event_number[(1],eventime({1])); */ 


/* find time delays by finding earliest channel event 
time, and subtracting that from the other channel 
times 

= 


firsttime 410000.0; 


for (l1=1;1<NCHAN;1++) { 
if(chan(1]) != 0 && event_flag{1l](1]) != 0 && 
chaneventtime(1])(1]) < firsttime) 
firsttime = chaneventtime(1) (1); 


} 


for (l=1;1<NCHAN; 1++) { 
if(chan[(1) != 0 && event flag[{1]{1]) != 0) { 
timedelay[(1] = ((chaneventtime(1)(1]) - firsttime) ; 
EPLInCe (pt eye teeta wets st Nh 4l,chan(1), 
timedelay{l], chaneventamp([1](1)); 


[= print to file to indicate end of event us 


for (l=1;l<num_active_events;1l++) { 
for (m=1; m<NCHAN;mt+) { 
event_flag(1](m] = event_flag[(1+1l] {m]; 
chaneventtime[1l][m]) = chaneventtime[1+1] [m)]; 
chaneventamp[(1])[m]) = chaneventamp(1+1)] (m]; 
} 
eventime(l]) = eventime(1l+1]; 
event_number(1l] = event_number([1l+1]; 
} 
num_active_events--; 
} 
time += 0.004*4.0; 
} 
} 
/* EXIT MODULE aL 


Seurce code for the Adetect proaram. 
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) = ecCOorineetSscderr, “processing record *#d\n", 
if(toggle == 1) { 
for(j=1;j3<NCHAN; j++) = { 
if(chan(3] != 0) { 
for (k=0;k<OBYTES; k++) 
oddobuf (3) (k] = 0; 


sq filt (sevenobuft[j] (0],soddobuft[j] (0],h, &tameseries([j] (0]); 


1_ave(j] = vconv(s&longfilt(0],&timeseries(j] (2* (FLN=-1)], 
LONGFILTLN) ; 


nskiptnseg); */ 


} 
} 
) 
else { 
for (j=1; J<NCHAN; j++) = [ 
re (chan [3], 0) { 
for (k=0;k<OBYTES; k++) 
evenobuf(j)(k]) = 0; 
sq_filt (&oddobuf(j] (0],&evenobuf(j] (0],h, &timeseries(j] (0)); 


1_ave[j) #*vconv(élongfilt [0], ctimeseries([j] (2*(FLN<-1)], 
LONGFILTLIN) ; 


Eor(k™0;K< (RECLN=2* (FLN=1));k+=4) { 
: for (i=l, NeHaN 1+) ei { 
1£ (chan[(1l] != 0) { 
Tiave(1] = (63.0*1 _ave[1l] + timeseries [1] (2* (FLN-1)+k]) /64.0; 
Woeaves=—s(Cimeserives (12> (= oN>1) +k) + 


timeseries([(1l] (2* (FLN-1)+k-1] + timeseries(1] (2* (FLN-1)+k-2] + 
timeseries(1] (2* (FLN-1)+k-3]) /4.0; 


/* reset old flags */ 


if (Elag(1) == ice 


(time - timeofflag(1]) > RESET_DELAY) 
flag (1) “<0; 


/* set flag if RATIO of signals is reached */ 


if ((sh_ave/1l_ave[(1])>=RATIO) [ 
1f (flag(ijleri1) { 


if (sh_ave > ampofflag(1l]) [ 
timeofflag(1l] = time; 
ampofflag{1] = sh_ave; 


} 


else { 
if (num_active_events == 0) [ 
flag({l] = 1; 
timeofflag(1] = time; 


Source cede for the Adetect program. 
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empont lag (i) = sn lave; 
} 
else { 
for (m=l;m<(num_active_events + 1);m+t+} { 
i= Kevent tlag(im] (lj -e=_1} { 
if ((time - chaneventtime[m] [1}) 
<= SIG DELAY) { 
if (sh_ave > chaneventamp[m] [1]) { 
chaneventtime(m](1] = time; 
chaneventamp([m] [1] = sh_ave; 
} 
} 
else { 
Elaq tly =. 13 
timeofflag(1] = time; 
ampofflag[(1l] = sh_ave; 
} 
} 


else { 


event ftlag({m] [1] =— 1; 
chaneventtime(m][({1] = time; 
chaneventamp[m] {1] = sh_ave; 
m = num_active_events; 


} 
/* end of set flag module */ 
/* start new event module */ 
flag_sum = 0; 


for (l=1;1<NCHAN; 1++) { 
if (flag(lj.S—.1) flag sum++; 

} 

if (((float) flag_sum/(float)nchan) >= THRESHOLD) { 
num_active_eventst+; 
eventime ([num_active_events] = time; 
number of events++; 
event_number(num_active_events] = number_of_events; 


for (1=1;1<NCHAN;1++) { 
event_flag(num_active_events][(1] = flag[1l]; 
chaneventtime (num_active_events]{1] = timeofflag[1]; 
chaneventamp[num_active_ events] {1} = ampofflag[1l]; 
flag({1l] = 0; 
timeofflag[{1] = 0; 
ampofflag(1] = 0; 
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} 
/* end of new event module */ 


/* start of deactivate old event module */ 


if (num_active events > 0 && (time - eventime({1]) > EVENT_DELAY) { 


forines (ptr, td ta tf\n", ZERO, 
event_number[(1),eventime(1]); 
l= fprintf (stderr,"td %f\n",event_number(1],eventime(1]); */ 


find time delays by finding earliest channel event 


/* 
time, and subtracting that from the other channel 
times 

a7 
firsttime — 10000.0; 
for (l=1;1<NCHAN;1++) { 

~ti(chan(1) != 0 4& event _flag{1]{1] != 0 && 
chaneventtime(1){(1] < firsttime) 

firsttime = chaneventtime(l1] [1]; 
} 
for (1=1;1<NCHAN;1++) = { 

LE(chan(1] != 0 && event _flag(1][1] != 0) { 
timedelay(1]) = ((chaneventtime(1l](1]) - firsttime); 
fprintr (per, “eC watcate sett \n',l,chan{l], 

timedelay(1l], chaneventamp([1] (1)); 
} 
} 
ees print to file to indicate end of event a7, 


for (l=1;l<num_active_events;1++) [{ 
for (m=l; m<NCHAN;m++) ({ 
event tlag(lj im] -= event flag({1+1] (m]; 
chaneventtime(1] ({m] = chaneventtime[(1+1] (m]; 
chaneventamp(1] (m] = chaneventamp[1+1)] [m]; 
} 
eventime(1l] = eventime(1+1]; 
eveneenumoce | Lie event number (1+1); 
} 


num_active_events--; 


} 


time += 0.004*4.0; 
} 
/* print out all events */ 


22 (numsactive events > 0) { 


Seurce code for the Adetect program. 
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for (k=1;k<(num_active_events + l);k++) { 
fprintf (ptr,"%d %d %f\n", ZERO, 
event_number([(k],eventime([k)); 


ig fprintf (stderr,"td %f\n",event_number[(k]),eventime(k]); */ 

/* find time delays by finding earliest channel event 
time, and subtracting that from the other channel 
times 

ey 


firsttime = 10000.0; 


for(l=1;1<NCHAN;1++) { 
if(chan(1) != 0 && event _flag[(k]){1l) != 0 && 
chaneventtime[(k) (1) < farsttime) 
firsttime = chaneventtime[(k] [1]; 


} 


for (1=1;1<NCHAN; 1++) { 
if(chan(1l) != 0 && event_flag(k] [1] != 0) { 


timedelay[{1l]) = ((chaneventtime([k]) [{1]) - firsttime) ; 


fprintf(ptr,"td %d $f %£\n",1,chan[l], 
timedelay{l], chaneventamp[(k) [(1)); 


} 
} 


* final summary to screen */ 


PpoEineLipts, td” ,--2Np) - 
Close (per) > 
fclose(iptr) ; 

exit(0) ; 


$q f£ilt (first, second, filter, output) 
/* filters and squares two data arrays */ 


£loat first [RECLN]), second (RECLN), f£ilter[FLN]; 
float output [RECLN+2*FLN] ; 


{ 
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PeepuL ecerOoS.in first (FLN =-1) places of output */ 


One 70 2s 
float transition(2*FLN]; 
float sum; 


for (i=0;i1<FLN-1; i++) 
Seutput [1] —. 0.0) > 


(= Puc an first: data 7/7 


for (L=FLN-1;i<RECLN;i++) { 
sum = yvoonv (&filter(0],éfirst(i}),FLN): 
output(i] = sqrt(sum * sum); 


} 


Pee oue in transition from first to second data */ 


for (jw0; 53<FLN-1;3++) { 
transition(j] = first (RECLN-(FLN-1) +3] ; 
transition[(j+FLN-1) = second[j]j; 


} 


for (i=0;1i<FLN-1; i++) { 

sum = veconv(éfilter(0),é&transition(i+FLN-1]) ,FLN) ; 
output (1+RECLN] = sqrt (sum * sum); 
) 


/*Soutein second data */ 
for (1=FLN-1;1<2*FLN; i++) { 
sum = vconv(&filter(0),&second(i),FLN); 


output (RECLN+i] = sqrt(sum * sum); 
) 
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Source Code for the farlecate Program 





HYSER‘'S GUIDE FOR THE location, farlocate, 
and finelocate PROGRAMS 

The purpose of the location programs is to find the 
spatial location of an event from the time delays between 
Signal arrival at different hydrophones. The program 
assumes a test location and computes the slant range to the 
individual hydrophones. The slant ranges are plotted 
against the experimental time delays and a least squares 
fit is done. The test location with the best least squares 
#1t 1S considered the location of the event. 

The input to the location programs is the output file 
of the detection program. Manual time delays may be 
substituted for the program generated time delays in this 
file, but this editting must be done before the location 
program is invoked. 

This location program is very user interactive. The 
user starts the location program with the command 

location “RETURN? 
The program asks for the input and output file names. Tt 
then reads the input file and asks whether the user would 
like to locate the first event. This allows the user to 
Skip down to the event of interest. The program then 
allows the user to adjust which hydrophone time delays will 
oe used in the location process. This is very handy for 
removing questionable time delays, in order to get a better 


Iscation solution. Wath the hydrophone channels chasen. 





oe tee 


the program proceeds with the actual location algorithm. 
In the finelicoate and fariecate programs the user is asked 
to specify which quadrant or direction is to be searched. 

The focation program tries test locations in a large 
Grid, and when the "best" location is found, then searches 
a smaller grid around this "best" location. The lecation 
and fTeriecate orograms have 4 levels of arids and the 
Ffinelacate oroaram has 3 levels. Intermediate answers are 
displayed for each level. 

TRe intermediate and final answers display the x and y 
coordinates of the best location, the standard deviation of 
the least squares fit (sigma), the group speed (which 
should be around 1440 m/s). and the y intercept of the time 
delay / slant ranqde plot. After the final answer the user 
1s asked whether or not he would like to remove outlying 
ooincts. If this option is selected the froaram removes 
Rycrosohoenes with a deviation from the least squares fit of 
more than 3 times the standard deviation, and the group 
speed and the standard deviation are recalculated and 
displayed. 

The user is then asked 1f he would like to locate the 
event with different Rydrophones, and if so returns the 
user to the start of the channel selection process. The 
“ser May repeat this ioccation scheme as many times az 
necessary for a particular event. Qnce the user 15 


ce: 


ial 


un 


fied with the lecation answer, and declines to locate 
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the event with different phones, the program calculates the 
event strancth based on the amplitudes in the input file, 
the location of the event and spherical spreading losses. 
The location parameters and the strength are then written 
to an outfile. 

At this point the user is given the option to exit the 
orogram or locate the next event in the input file. Tre 
location process continues until the user exits or until 


the last event ina fileis iocated. 





eo 


/* location.c 


program to locate the spatial position of an event based on 
time delays taken from “detection”. 


Source strength is also computed. 


mi 


#include <stdio.h> 
#include <math.h> 


#define 
#define 
#define 
#define 
#define 


main () 
{ 


floae 
float 
Eloat 
float 
£loat 
float 
float 


ime, 


PHONES 31 

LEVEL 4 

FINENESS 20 

DEPTH 91.0 

Satioliayiay 0; 0000000112202 


amo [PHONES], phonex(PHONES], phoney [PHONES]; 
timedelay (PHONES), r{PHONES], bestrange [PHONES]; 
sumtime, sumr, sumrsq, sumtimer, slope, yintr; 

Sigma, bestsigma, bestslope, bestyintr; 

bestamp, N, source(PHONES), sumsource; 

as, yY9S, XCMtr, ycontr, a, b, bestx, besty:; 

level, gridsize, xfineness, yfineness, time, gpspeed; 


j, tape, date, hour, min, sec, event, flag; 


int phone, num, 1, n, m, bestflag; 
int phoneflag(PHONES]; 
int eventselect, rerun, change, answer, bye, quadrant; 


char iname(80], oname([(80]; 


PLlLes pcr, Open(); “Llocptr, *optr; 


/* PROGRAM INITIALIZATION */ 


bye = 


OF 


fprineedstderr,"input file = \n"); 
scanf("%s", iname) ; 
fpeints (tderr, "output file = \n"™); 
scanf("%s", oname); 


/* open 


files */ 


if((ptr = fopen(iname, “r")) == NULL) { 
fprinca (stderr, “can’t open ts\n", iname) ; 
Cxit (1); 


Source code for the farlocate program 
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if((optr = fopen(oname, "w")) == NULL) { 
fprintf(stderr, "can’t open %s\n", oname) ; 
exit (1) ; 

} 


if((locptr = fopen("array_loc", "r")) == NULL)  { 
fprintf(stderr,"can’t open array_loc file\n"); 
exit (1); 

} : 


/* read hydrophone locations into array */ 


for (i=1; i<PHONES; i++) ( 
fscanf(locptr,"%d %f %f£",&phone, &phonex[i], &phoney([i)); 


} 


/* read input file header */ 


fscanf(ptr, "td td td %d %d", &tape, f&date, G&hour, &min, &sec); 
/* read event header */ 


fscant (ptr, "3d ", “eflacg); 
while(bye != -1)  { 


eventselect = 0; 
while(eventselect '= 1) { 


Lttftag <0) 
exit (0); 


for (i=; i<PHONES; i++) { 
phoneflag(i} = 0; 
} 


fscanf(ptr, “sd %£", &event, &time); 


for(i=l; i<PHONES; i++) { 
£scant (ptr, "ta", &flaq) ; 


PeCelag > 0). { 
Esecant (ptr, “=*d", &5); 
phonetilag(j )*= 1; 
fscanf(ptr, “tf %£", &timedelay(j]), samp(j)});: 
} 
else { 
i = PHONES; 
} 
} 


fprintf(stderr, "event = %d, time = %f\n", event, time); 


Source code for the ferlocate program 
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fEprintcsederr, Do you wish to locate this event? (1 = yes) \n"); 
scanf("%d", &eventselect); 
} 


/* channel selection */ 


rerun = 1; 

while (rerun == 1) [ 
change = 1; 
while (change == 1) /[ 


perme. (Sstders, phone delay Nac) 
for(i=l; i<PHONES; i++) ( 
irf(phoneflag({i] '= 0) ( 
fprintt (stderr, ™ td $£ \n",i,timedelay([i]J); 


} 
} 


fprinté (stderr,”"Do you wish to change status? (1 = yes)\n"); 
scanf ("%d", é&change) ; 
if (change == 1) { 


fprintf(stderr, "change status by typing phone#\n") ; 
fprintf(stderr, "type -1 to quit\n"™); 
scant ("%td",&3); 


while (3 != -1) ([{ ‘ 

L£ (phoneflag(j3] != 0) { 
phoneflag{j] = 0; 

} 

else [{ 
phoneflag{j] = 1; 

} 

seanrt. sad", &3); 


} 


/* locate event */ 


num = QO; 
sumtime = 0.0; 


for (imsl; 2<PHONES; i++) ( 
if (phoneflag({i] != 0) ({ 
numt+; 
sumtime += timedelay(i]; 
} 
} 
fprintf (stderr,”"select quadrant to search (1=NE, 2=NW, 3=SW, 4=SE\n"); 
fprintf(stderr, ”"5=N, 6=S, 7=E, 8=W)\n"); 
scanf ("td", &quadrant); 
if (quadrant == 1) [ 


Source code for the farlocate program 





xenter 
Vencs 

} 

else if 
xener 
Yontc 

} 

else if 
xeGntr 
Vener 

} 

else if 
Senter 
Vents 

} 

else if 
xentc 
Vent © 

} 

else if 
xecntxc 
Ventr 

} 

else if 
xente 
Vener 

} 

else if 
Sent £ 
Venter 


} 


else { 
xentr 
yentr 


} 


= 10000.0; 
=~ 10000 .0; 


(quadrant == 
= ~10000.0; 
= 10000.0; 


(quadrant == 
= ~10000.0; 
= ~10000.0; 


(quadrant == 
= 10000.0; 
= -10000.0; 


(quadrant == 
=. 0.0; 
= 10000.0; 


(quadrant == 
= 0.0; 
a =1 000020; 


(quadrant == 
= 10000.0; 
= 0.0; 


(quadrant == 
= ~10000.0; 
= 0.07 


=e 


bestflag = 0; 


2) 


3) 


4) 


3) 


6) 


7) 


8) 


ee 


for(1l=0; 1<LEVEL; l++) { 
level = LEVEL-1-1; 
gridsize = pow(10.0,level); 


for(m=0; m<FINENESS; m++) { 
yfineness =m ~- (FINENESS/2); 
ygs = yentr + yfineness * gridsize; 


for(n=0; n<FINENESS; nt++) ( 
xfineness = n ~- (FINENESS/2); 
xgs = xcntr + xfineness * gridsize; 
sumr = 0.0; 
sumzsq = 0.0; 
sumtimer = 0.0; 
for (i=l; i<PHONES; i++) { 
if(phoneflag{i) !'= 0) { 
a = xgs~-phonex[i]; 


Source cede for the farlocate Droaqram 
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b = ygs-phoney (iJ; 
Elz) = saqcst (pow(a,2.0) + pow(b, 2.0) + pow(DEPTH, 2.0)); 


sumx += r{i); 
sumrsq += pow(r[(1],2.0); 
sumtimer += (r{ijJ * timedelay({(i)); 


N = num; 

slope = ((N * sumtimer) - (sumr * sumtime) ) / 
CUCIN  *Ssumreq) = pow (sumr, 2.0) ) ; 

yintr = (sumtime - (slope * sumr))/N; 


sigma = 0.0; 
for(isml; i<PHONES; i++) { 
if(phoneflag{iJ != 0) { 
sigma += pow((timedelay(ij~-yintr-(slope * rf{iJ)),2.0); 
} 
} 
sigma = sqrt(sigma/N); 


if (bestflag == 0) { 
bestx = xg3s; 
besty = ygs; 
bestsigma = sigma; 
bestslope = slope; 
bescvintr = vintx©; 
for (il; i<PHONES; i++) § { 
1£f (phoneflag != 0) { 
bestrange[{ij] = r{ij; 
} 
} 
bestflag = 1; 
} 


else { 
if (sigma < bestsigma) § { 
bestx = xgs; 
besty = ygs; 
bestsigma = sigma; 
bestslope = slope; 
bestyintr = yintr; 
for (i=l; 1i<PHONES; i++) § { 
Lf “fohenef lag '=.0) ~ { 
bestrange({ij] = r{ij; 


seurce code for the farlocate program 
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gpspeed = 1.0/bestslope; 


forinee (stderr, "bestx = %£, besty = %£, sigma = %£\n", 
bestx, besty, bestsigma); 

fprints (stderr,"group velocity = %tf\n", gpspeed); 

Eprintc (Scdcrr, y intercept = tf\n", bestyintr) ; 

xcntr = bestx; 

yentrF= besty:; 


fprintf (stderr,"Do you wish to remove outlying points? (l=yes)\n"); 
scan£t ("td", answer) ; 
if (answer == 1) [{ 

num = 0; 

sumr =0.; 

sumtime = 0.; 

sumrsq = 0.; 

sumtimer = 0.; 


for (i=l; i<PHONES; i++) { 
L£{pnonetiag(2) '— 0) { 
if (sqrt (pow((timedelay [i] -bestyintr- (bestslope*bestrange[i]})), 
2.0)) < 2.S5*bestsigma) { 
numn++; 
sumtime += timedelay([(i]; 
sumr += bestrange([i]; 
sumrsq += bestrange([(i]) *bestrange [i]; 
sumtimer += bestrange(i] *timedelay[i]); 
} 
else [ 
forint=|(Stcer=r, CutlVing onone .*..sa\n", i); 
} 
} 
} 
N = num; 
bestslope = ((N*sumtimer) - (sumr*sumtime) ) / 
CU@='sunz sq) ~ (Sums *sumz) ) : 
bestyintr = (sumtime-(bestslope*sumr) )/N; 
gpspeed = 1.0/bestslope; 
fnorint. (Stderr, “bestx = %£, besty = %£,. Sigma ™ %£\n", 
bestx, besty, bestsigma); 
Sprints (stderr, group velocity .= 4£\n", gpspeed); 
ferinee (stderr, y intercept, “sttin , bestyintr) ; 
} 


fprintt (stderr, 


"Do you wish to relocate with different phones? (l=yes)\n"); 
Ream t sa ;& rerun) ; 


source code for the farflocate program 
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/* €inding source amplitude */ 


for(i=l; i<PHONES; i++) { 
if (phoneflag[i] != 0) { 
source[i] = (amp({i]}/SENSITIVITY) * bestrange[iJ; 
sumsource += source[i)> 
} 
} 
bestamp = sumsource/N; 


ESrince (optr, “td td S£& S£ FE SE\n", event, date, time, 
bestx, besty, bestamp); 
fprintf(optr,"t£ t£ t£\n", bestsigma, gpspeed, bestyintr); 


Eprintel (stderr, "Do you wish to exit? (-1 = ex2t)\n"); 
scant ("%d", é&bye); 
} 
SClLOsSe4{ptr) ; 
EGLGSe(Optr) ; 
Eclose(locptr); 
ex2e (0); 


} 


Source code for the farlocate program 
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APPENDIX C 


rae ©—1: Event Location Summary 


Table C-2: Tape Summary 


Table C-2: Event Interarrival Time Summary 


Table C-4: Event Strength Summary 
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Tacle C-1 Event Location Summary 
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Waae ¥ ! f R phi 51QMa sigma c 0 Ry Graprone: 
(jane F, Reraved 
event #) (ai) ia! (a) (m) idecrees) (sec) (fi) (m/sec) 
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Table 0-2 Spreading Loss Function, Gir) 


Theta zero Theta | z Ri R2 G(R1) G{R2) 


10 log = 16 log 
G(RL) G(R2) 


Theta zero is the surface launch angle. 

Theta 1 15 the angle of the refractive path at the hydrophone depth (93m) 

z 1s the aaxiaue depth of the refractive path. 

R1 is the horizontal range to the hydrophone intersected on the downward swing 
R2 15 the horizontal range to the hydrophone intersected on the upward swing 
G(R) 1s the spreading loss function. 
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0.085 0.916 96 1839 2697 4.2€+08 6.2E+08 86.239 87.901 
0.086 0.021 99 1741 2848 3.9E+08 6.3E+08 85.881 88.018 
).087 0.029 101 1666 2977 3.6E+08 6.9608 83.975 88.097 
0). 088 0.028 103 1603 3094 3.4E+08 6.SE+08 85.502 88.197 
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